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Abstract 

This  report  outlines  the  initial  findings  of  the  research  team  conducting 
the  ERDC-directed  research  project  “Nanoscale  Studies  of  Polycrystalline 
Materials  with  Emphasis  on  Ceramics  Syntheses.”  It  provides  an  assess¬ 
ment  of  the  state-of-the-art  in  the  multi-scale  simulation  methods  that  can 
predict  polycrystalline  ceramic  mechanical  properties  and  ceramic  sinter¬ 
ing  from  basic  physics  and  material  structure.  The  report’s  findings  will  be 
used  to  identify  strengths  and  weaknesses  in  the  technology,  to  under¬ 
stand  how  the  different  simulation  components  must  fit  together,  and  to 
guide  follow-on  research  programs  towards  the  long-term  development  a 
ceramic  composite  that  has  fracture  toughness  and  tensile  strength  ap¬ 
proximately  5  times  that  of  existing  polycrystalline  ceramics,  such  as  sili¬ 
con  carbide  or  boron  carbide.  If  such  a  ceramic  composite  were  developed, 
then,  based  on  current  strength-to-weight  and  stiffness-to-weight  ratios,  it 
could  replace  steel  and  aluminum  for  most  structural  applications  with  an 
attendant  two-thirds  reduction  in  weight.  This  would  have  enormous  im¬ 
pact  on  Army  portable  protective  structures,  equipment,  and  logistics.  Key 
to  this  development  is  the  growing  capability  in  numerical  simulations  to 
predict  material  behavior  based  on  atomic  and  crystalline  morphology. 
Such  simulations  provide  new  insight  into  the  causal  relationships  be¬ 
tween  material  structure  and  material  behavior.  The  simulations  can 
guide  both  polycrystalline  material  design  and  synthesis  methods  such  as 
sintering. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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Preface 

This  report  provides  an  assessment  of  the  state-of-the  art  in  multi-scale 
numerical  simulation  methods  necessary  to  model  the  mechanical 
properties  and  the  synthesis  of  polycrystalline  ceramics.  It  identifies 
strengths  and  weaknesses  in  these  technologies  and  how  the  methods  can 
be  used  synergistically  to  develop  a  fundamental  understanding  of  the 
underlying  phenomena.  The  findings  of  this  report  will  be  used  to  guide 
ERDC’s  materials  development  research  whose  long-term  goal  is  the 
development  of  a  structural  ceramic  that  can  replace  steel  and  aluminum 
for  most  applications  with  significant  (two-thirds)  weight  reduction. 

This  report  is  a  product  of  the  ERDC-directed  research  program 
“Nanoscale  Studies  of  Polycrystalline  Materials  with  Emphasis  on 
Ceramics  Syntheses.”  This  program  is  a  joint  ITL-GSL-CERL-EL  program 
under  co-leads  Dr.  Charles  R.  Welch  (ITL)  and  Dr.  John  F.  Peters  (GSL). 
The  report  was  prepared  by  the  research  team  members  identified  and  was 
edited  by  Dr.  John  Peters.  Marsha  Gay  and  Holly  Clavien  provided 
technical  editing  and  formatting.  The  report  was  prepared  under  the 
general  direction  of  the  Laboratory  Directors  of  the  four  ERDC 
laboratories:  Dr.  Reed  Mosher  (ITL),  Dr.  David  Pittman  (GSL),  Dr.  Ilker 
Adiguzel  (CERL),  and  Dr.  Beth  Flemming  (EL),  and  under  the  general 
supervision  of  the  GSL  Technical  Director,  Dr.  David  Horner  who  oversaw 
ERDC’s  material  development  efforts. 

COL  Kevin  Wilson  was  the  Commander  of  ERDC  during  the  time  of  this 
report;  the  Associate  Director  of  ERDC  was  Dr.  John  Cullinane;  and  the 
ERDC  Executive  Director  was  Dr.  Jeffery  P.  Holland. 
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1  Introduction 

This  report  outlines  the  findings  of  the  initial  phase  of  the  project 
“Nanoscale  Studies  of  Polycrystalline  Materials  with  Emphasis  on 
Ceramics  Syntheses.”  The  objective  of  the  project  is  to  develop  multiscale 
modeling  techniques  for  polycrystalline  ceramic  sintering,  principally  the 
sintering  of  silicon  carbide  (SiC).  This  research  advances  the  long-term 
effort  to  increase  the  tensile  strength  and  fracture  toughness  of 
polycrystalline  ceramics  to  a  factor  of  approximately  five  that  of  existing 
values.  Sintering  is  the  process  by  which  polycrystalline  ceramics  are  made 
from  crystalline  powder  by  applying  heat,  stress,  and  sometimes  electric 
and  magnetic  fields.  Polycrystalline  ceramics  offer  the  promise  of  superior 
structural  materials  having  strength-to-weight  and  stiffness-to-weight 
ratios  that  are  factors  of  three  or  more  those  of  current  structural 
materials  such  as  aluminum  and  steel.  In  addition,  ceramics  retain  their 
properties  under  high  temperature  (over  3000  °F),1  resist  corrosion,  and 
are  created  from  abundant  materials.  Unfortunately,  bulk  specimens  of 
ceramic  materials  have  low  fracture  toughness  and  tensile  strength 
compared  with  the  structural  metals,  making  them  unsuitable  for  many 
important  structural  applications.  Moreover,  it  has  been  the  experience  of 
the  research  community  that  the  mechanical  properties  of  ceramics  have 
not  been  improved  significantly  using  traditional  empirical  approaches  of 
manipulating  composition,  particle  characteristics,  and  sintering 
properties.  A  deeper  understanding  of  fundamental  processes  is  required. 

Technical  ceramics  have  been  around  for  over  a  century.  SiC,  for  example, 
was  manufactured  and  sold  by  Edward  Acheson  as  early  as  1893.  It  can  be 
made  by  combining  silicon  sand  and  carbon  in  a  graphite  electric  resis¬ 
tance  furnace  at  temperatures  of  1600  to  2500  °C.  SiC  has  a  density  of 
194  lb/ft3,  compressive  strength  of  570  ksi,  Young’s  modulus  of 
59.5  million  psi,  a  flexure  strength  (related  to  tensile  strength)  of  80  ksi, 
and  a  fracture  toughness  of  4.6  MPa  m°-5. 

Fracture  toughness  is  related  to  brittle  behavior— the  lower  the  fracture 
toughness,  the  more  brittle  the  material.  For  comparison,  high-strength 
7075-T6  aluminum  has  a  fracture  toughness  of  24  MPa  m°-5;  4340  steel,  a 


1 A  table  of  factors  for  converting  non_SI  units  of  measurement  to  SI  units  can  be  found  on  page  ix. 
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very  high  strength,  very  tough  steel,  has  a  fracture  toughness  of  50  MPa 
m°-5;  while  concrete  has  a  fracture  toughness  of  0.2  to  1.4  MPa  m°-5. 


Table  1  compares  several  common  ceramics  with  high-strength  aluminum 
and  very  high-strength  steel.  The  ceramics  have  densities  roughly  that  of 
aluminum;  Young’s  moduli  from  four  to  six  times  that  of  aluminum,  and 
compressive  strengths  from  four  to  seven  times  that  of  aluminum.  The 
densities  of  the  ceramics  are  from  one-third  to  one-half  that  of  the  steel; 
Young’s  modulus  from  1.5  to  2  times  that  of  steel;  and  compressive 
strength  1.5  to  3  times  that  of  very  high  strength  4340  steel  alloy.  The 
extreme  values  for  ceramics  are  the  reasons  they  are  used  for  tool  cutting, 
for  high-temperature  applications,  for  aerospace  applications,  for  bearing 
surfaces,  for  precision  parts,  and  for  armor,  among  other  things.  Their  low 
tensile  strength  and  low  fracture  toughness  prevent  their  use  in  structures, 
bridges,  vehicles,  airplanes,  ships,  trains,  gear  sets,  turbines,  etc. 


Table  1.  Material  properties  of  ceramics  compared  with  steel  and  aluminum. 


Material 

Young’s 

Modulus 

(Million 

psi) 

Compressive 

Strength 

(ksi) 

Density 

(lb/ft3) 

Fracture 
Toughness 
(MPa  m0-5) 

Bending 

Strength 

(ksi) 

Maximum 
Service 
Temperature 
(F  Degrees) 

Silicon 

Carbide 

(SiC) 

59.5 

566 

194 

4.6 

80 

3000 

Boron 

Carbide 

(B4C) 

63.8 

560 

157 

3.1 

62 

1112 

Aluminum 

Oxide 

(AI2O3) 

43.5 

305 

230 

3.5 

47 

3092 

High- 

Strength 

Steel 

(4340) 

29.7 

160  -  200 

490 

50 

160  to 

200 

(Melting  Point 
-  2600) 

Aluminum 

7075-T6 

10.4 

74-78 

170 

24 

74-78 

(Melting  Point 
-  900) 

It  is  hypothesized  (Welch  et  al.,  2009)  that  if  a  ceramic  matrix  is 
appropriately  married  to  carbon  nanotubes  or  to  graphene,  the  carbon 
nanotubes  or  graphene  would  provide  the  tensile  strength  and  tensile 
load-carrying  capability  to  overcome  the  fracture  toughness  and  tensile 
weaknesses  in  ceramics.  It  is  also  hypothesized  that  such  a  composite 
would  operate  on  a  nanoscale  much  like  steel-reinforced  concrete  operates 
on  a  macroscale,  arresting  cracking  before  it  can  propagate. 
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The  technical  literature  provides  encouraging  results  to  support  this 
hypothesis.  The  work  for  the  most  part  is  empirical.  Karandikar  et  al. 
(2008)  incorporated  carbon  fiber  and  carbon  nanotubes  into  SiC  and 
boron  carbide  matrices.  They  report  an  improvement  in  the  SiC  fracture 
toughness  of  from  4  to  7  MPa  m°-5,  a  73  percent  increase,  although  still 
below  the  long-term  goals  of  the  program.  Yamamoto  et  al.  (2008)  built  a 
highly  homogenous  composite  of  alumina  and  multiwalled  carbon 
nanotubes  that  was  found  to  exhibit  a  27  percent  improvement  in  bending 
strength  and  a  25  percent  improvement  in  fracture  toughness  to  100  ksi 
and  5.0  MPa  m°-5,  respectively.  From  the  work  of  Steif  and  Trojnacki 
(1992)  on  ceramics,  it  can  be  inferred  that  increases  in  bending  stress 
correspond  to  improvements  in  postyield  behavior.  Hence  both  parameter 
improvements  reported  by  Yamamoto  et  al.  suggest  reductions  in  brittle 
behavior  of  the  alumina  as  a  result  of  the  multiwalled  carbon  nanotubes. 
Kwon  et  al.  (2009)  reported  aluminum  (Al) /carbon  nanotube  (CNT) 
composites  with  tensile  strengths  that  were  three  times  larger  than  pure 
aluminum  because  of  the  CNT  reinforcement. 

Xia  et  al.  (2004)  constructed  a  highly  ordered  multiwalled  CNT  Al 
composite  and  subjected  it  to  nanoindentation.  Their  analyses  determined 
that  the  multiwalled  CNT  members  improved  the  fracture  toughness  of  the 
material  by  5  Mpa-m°-5.  They  also  analyzed  the  apparent  tensile  strength 
of  the  CNTs  in  the  matrix  and  the  interfacial  stresses  between  matrix  and 
the  CNTs  and  determined  values  of  15-25  GPa  and  40-200  MPa, 
respectively,  for  these  parameters.  Wang  et  al.  (2006)  fabricated  a  CNT- 
SiC  composite  using  a  high-pressure  (2  to  8  GPa),  high-temperature 
(1770-1970  K)  reactive  sintering  process.  They  report  a  fracture  toughness 
of  6. 8-7.1  MPa-m°-5for  the  composite,  well  above  the  nominal  values  for 
SiC. 

CNTs  have  improved  mechanical  performance  in  the  cases  cited.  However, 
in  some  cases  the  addition  of  CNTs  has  produced  marginal  or  inferior 
mechanical  properties,  as  discussed  in  Zhan  et  al.  (2003).  One  potential 
problem  with  adding  CNTs  is  the  difficulty  in  dispersing  them  completely 
throughout  the  SiC  matrix.  Without  dispersal,  the  CNTs  clump  as  a  result 
of  van  der  Waals  forces  between  the  molecules,  ultimately  leading  to  voids 
in  the  final  polycrystalline  matrix  that  become  sites  where  fractures 
initiate  to  cause  brittle  failure.  Evaluating  reports  of  the  success  or  failure 
of  adding  CNTs  to  enhance  strength  and  toughness  of  ceramic  composites 
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requires  knowing  details  of  how  the  CNTs  were  incorporated  into  the 
composite  and  whether  they  were  functionalized  beforehand. 

The  reported  gains  in  mechanical  properties  indicate  the  following: 
performance  improvements  can  be  achieved  with  CNT-enhanced 
composites,  the  mechanisms  contributing  to  these  improvements  are 
unclear,  and  the  fundamental  principles  to  guide  improvements  of  the 
magnitude  required  to  meet  the  longterm  goals  of  this  research  are 
lacking.  Ceramic  researchers  have  recognized  the  twin  weaknesses  of  the 
low  tensile  strength  and  brittle  behavior  of  ceramics  for  decades.  Their 
approach  to  improving  ceramics,  as  in  the  studies  cited,  has  been  largely 
empirical,  with  the  development  being  hit-and-miss  and  consisting  of 
cycles  of  build  then  test;  rebuild,  test,  etc.  Their  testing  approach  has 
largely  been  at  the  macroscale.  This  development  approach  is  common  to 
the  historical  development  taken  by  most  other  material  developers. 
Material  response  usually  begins  at  the  nanometer  scale.  Using  only 
empirical  approaches  to  material  response  often  precludes  the 
development  of  important  understanding  of  underlying  phenomena. 
Relying  on  macroscale  experiments  exacerbates  the  problem  as  details  of 
important  nanoscale  phenomena  remain  hidden. 

The  use  of  atomistic  simulations  to  understand  material  behavior  at  the 
nanoscale,  coupled  with  coarse-grained  simulations  to  propagate  that 
understanding  across  larger  material  volumes,  provides  critical  insight 
into  what  causes  a  particular  material  response.  Such  simulations  must  be 
validated— and  in  some  cases  calibrated— with  experimental  evidence 
across  the  same  scales.  Experiments  likewise  have  error,  either  systematic 
or  random,  and  leave  gaps  in  defining  important  quantities.  Simulations 
help  to  understand  the  experiments,  define  physical  quantities  that  cannot 
be  measured,  and  help  identify  experimental  error.  Thus  experiments  and 
simulations  complement  one  another.  Both  nanoscale  material 
diagnostics  and  response  and  multiscale  material  simulations  are 
developing  technologies  that  are  being  used  as  they  are  being  developed. 
Used  together  they  can  change  material  development  from  being  an 
evolutionary  process  to  a  revolutionary  process. 

1.1  Multiscale  design 

Heterogeneous  composites  such  as  polycrystalline  ceramics  must  be 
understood  from  the  perspective  of  various  scales.  The  fundamental 
properties  of  the  crystalline  materials  that  make  up  the  composites  are 
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determined  at  the  atomic  scale.  However,  the  engineering  properties  are 
determined  as  much  by  the  interaction  among  crystalline  grains  as  the 
strength  of  the  grains  themselves.  Thus,  a  multiscale  approach  must  be 
considered.  In  this  report,  five  scales  will  be  considered,  based  on  the 
methods  used  to  investigate  the  physics: 

•  At  the  smallest  scale,  quantum  mechanics  is  needed  to  understand 
atomic  interactions. 

•  At  the  molecular  scale,  Newtonian  mechanics  is  applicable  using  po¬ 
tential  functions  derived  from  quantum-scale  studies  to  define  intera¬ 
tomic  forces  (e.g.,  molecular  dynamics). 

•  Although  a  great  deal  can  be  learned  about  structural  response  using 
molecular  dynamics  (MD),  the  time  scale  of  mechanical  interactions  is 
much  shorter  than  that  of  diffusive  mechanisms,  which  are  important 
in  considering  material  synthesis.  For  many  processes,  MD  is  simply 
not  feasible  and  Monte  Carlo  (MC)  methods,  which  determine  equilib¬ 
rium  positions  of  atoms,  must  be  used. 

•  At  length  scales  of  greater  than  a  few  nanometers,  atomic  models  give 
way  to  mesoscale  models  that  depict  the  material  as  discrete  entities 
whose  interactions  are  described  through  presumptive  laws,  either  de¬ 
rived  from  MD  computations  or  directly  observed  in  nanoscale  experi¬ 
ments. 

•  The  mesoscale  extends  into  the  micron  range  where  so-called  continu¬ 
um  scale  descriptions  are  used.  The  continuum  scale  refers  to  the  abil¬ 
ity  to  define  a  representative  volume  element  that  is  large  enough  that 
mechanical  interactions  between  discrete  elements  of  the  mesoscale 
can  be  homogenized  into  equivalent  continuum  properties,  which  then 
can  be  endowed  to  infinitesimal  elements  for  purposes  of  engineering 
analysis. 

While  not  formally  considered  here,  there  is  also  a  structural  scale  in 
which  the  material  can  be  idealized  as  structural  elements  (plates,  beams, 
etc.)  whose  mechanics  are  derived  from  approximations  to  the  continuum 
equations. 

This  report  approaches  the  multiscale  problem  as  consisting  of  the 
following  steps: 

•  Given  target  properties,  a  mesoscale  model  is  developed  of  the  grain 
structure.  This  model  spans  the  range  of  10  nm-i  pm.  In  this  step  the 
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optimal  grain  structure  and  required  grain  boundary  properties  are 
determined. 

•  Molecular  dynamics  models  of  the  grain  boundary  are  developed  to 
determine  the  structure  required  to  achieve  the  properties  required. 
The  scales  involved  are  be  less  than  too  nm. 

•  The  feasible  chemical  structures  that  provide  suitable  mechanical 
properties  of  grain  boundaries  are  determined  from  Monte  Carlo 
simulations  and  ab  initio  computations.  These  simulations  will  define 
the  pressure-temperature  states  required  to  synthesize  the  materials. 

•  Continuum-scale  models  are  used  to  examine  the  temperature- 
pressure  regime  within  the  sintering  chamber.  The  scales  of  these 
models  span  the  dimensions  of  the  sintering  device. 

•  The  characteristics  of  the  green  body  (presintered)  required  to  achieve 
the  final  material  structure  are  evaluated  using  experiments  and 
methods  similar  to  granular  media.  The  scale  for  these  investigations  is 
similar  to  that  of  the  mesoscale  modeling. 

This  process  will  be  modified  as  necessary  to  accommodate  developments 
in  the  research.  Note  that  the  process  is  neither  a  top-down  nor  a  bottom- 
up  sequence  as  is  often  depicted  in  descriptions  of  multiscale  analysis. 
Rather,  it  reflects  the  process  of  material  design  in  which  material 
response  must  be  considered  in  tandem  with  material  synthesis.  The 
properties  of  the  homogenized  medium  are  derived  from  interactions  at 
the  mesoscale.  New  concepts  of  nanostructure,  derived  from  studies  of 
biological  materials  (Buehler  2010;  Garcia  et  al.  2010;  Bratzel  et  al.  2010; 
Espinosa  et  al.  2011),  provide  avenues  for  advancement  in  performance 
that  have  been  exploited  but  to  a  limited  extent  in  ceramics  research  (Rim 
et  al.  2011).  Mesoscale  analyses  can  be  used  to  determine  the  mechanical 
behavior  of  proposed  structures.  Once  the  pressure  and  temperature 
environment  favorable  to  the  structure  is  determined,  the  issue  of  creating 
those  conditions  at  a  practical  size  must  be  addressed  through  continuum- 
scale  computations.  Finally,  issues  such  as  how  the  green  body  of  granular 
powder  is  to  be  created  and  how  it  is  to  be  modified  through  sintering 
must  be  addressed  through  mesoscale  models.  Beyond  this  modeling 
process,  many  practical  considerations  of  material  synthesis  and 
composition  can  be  addressed  only  through  laboratory  experimentation 
combined  with  material  diagnostics. 
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1.2  Outline  of  report 

The  report  is  presented  in  six  chapters  in  the  approximate  order  of  the  five 
steps  discussed  in  the  previous  section.  The  report  is  the  product  of  the 
initial  studies  with  the  expectation  that  the  report  will  be  revised  as  need¬ 
ed  as  the  technology  and  research  progresses.  The  report  is  not  a  detailed 
description  of  the  methods  presented;  rather  documents  the  methods  in¬ 
vestigated,  technical  gaps  identified,  and  the  rationale  for  future  research 
tasks. 

1.3  Conclusions  of  2010  effort 

The  technology  necessary  for  implementing  the  five-step  procedure  for 
modeling  improved  SiC-based  polycrystalline  materials  has  been  demon¬ 
strated,  although  several  technical  gaps  remain: 

•  At  the  mesoscale  the  mechanical  response  of  polycrystalline  materials 
can  be  modeled  up  to  the  point  of  failure.  Fracture  along  grain  bound¬ 
aries  can  be  modeled  with  finite  element  technology  although 

trans  granular  failure  is  more  difficult.  Methods  proposed  for  modeling 
fracture,  such  as  peridynamics,  were  investigated  but  their  underlying 
physical  foundations  were  unclear.  The  present  capability  is  sufficient 
to  compare  polycrystalline  assemblages  in  terms  of  resistance  to 
intergranular  failure. 

•  Molecular  dynamics  is  applicable  to  assessing  the  structural  response 
across  grain  boundaries  and  the  effects  of  inserting  CNT  and  graphene 
within  the  SiC  matrix.  MD  is  not  suitable  for  modeling  the  synthesis 
process  per  se  because  of  extreme  mismatch  between  the  time-steps 
that  must  be  used  in  a  stable  MD  computation  and  the  computational 
duration  required  for  the  synthesis  processes  to  play  out.  A  critical 
element  of  the  MD  work  is  to  find  potential  functions  appropriate  for 
each  structure  investigated. 

•  The  sintering  process  can  be  captured  with  the  so-called  Potts  model, 
although  much  of  this  work  in  the  literature  has  been  done  in  two 
dimensions.  To  be  useful  for  the  present  work,  three-dimensional 
(3-D)  capabilities  are  needed.  Monte  Carlo  methods  can  possibly 
provide  critical  information  and  will  be  the  subject  of  further 
investigations. 

•  Once  the  ideal  sintering  conditions  are  determined,  the  ability  to 
obtain  these  conditions  in  full-scale  sintering  devices  can  be  evaluated 
using  continuum  modeling.  A  preliminary  multiphysics  analysis  was 
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performed  to  demonstrate  this  capability,  which  was  documented  in  a 
conference  presentation  and  paper  (Allen  et  al.,  2011). 

•  A  number  of  sintering  methods  were  reviewed.  The  Spark  Plasma 
Sintering  (SPS) ,  which  is  also  known  as  the  Field  Assisted  Sintering 
Technology  (FAST)  was  identified  as  the  focus  of  future  laboratory 
efforts  because  SPS  with  its  attendant  lower  sintering  temperatures 
and  faster  sintering  times  was  successful  in  producing  SiC-CNT 
composites  in  which  the  CNTs  remained  intact  (Inam  et  al.,  2010).  One 
intriguing  aspect  of  SPS  and  CNTs  is  the  apparent  ability  of  the  CNTs 
to  act  as  a  sintering  aid,  producing  smaller  crystalline  grains  in  the 
composite. 

An  impression  gained  from  the  work  of  2010  is  that  beyond  a  brute-force 
reinforcing  of  grain  boundaries  using  CNT  or  graphene,  the  goal  of  in¬ 
creasing  both  strength  and  toughness  can  be  met  only  by  manipulating  the 
polycrystalline  structure.  Previous  studies  using  particle  methods  to  study 
the  nanogranular  behavior  of  C-S-H  (C=CaO,  H=H20,  S=Si02)  in  ce¬ 
ment  showed  the  brittle  nature  of  equidimensional  (nearly  spherical) 
grains  (Chandler  et  al.  2009).  Comparison  to  biological  materials  further 
suggests  that  composites  of  strong  and  weak  components  are  required  for 
toughness.  Thus,  key  questions  are  what  structure  is  suitable  and  how 
might  that  structure  be  obtained? 
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2  Atomic-Scale  Modeling 

The  purpose  of  atomic-scale  modeling  in  this  research  is  to  use  atomistic- 
scale  and  mesoscale  simulation  to  study  materials  processing  by  linking 
the  performance  of  high-performance  materials  with  the  thermo¬ 
mechanical  processing  steps  required  to  produce  them.  This  chapter  dis¬ 
cusses  the  atomistic  modeling  and  simulation  progress  of  polycrystalline 
materials  synthesis  by  researchers  at  the  U.S.  Army  Engineer  Research 
and  Development  Center  to  date,  which  has  focused  on  identifying  the  re¬ 
quirements  for  modeling  the  synthesis  of  ceramic  polycrystalline  materials 
at  the  atomic  level.  This  chapter  also  presents  plans  for  future  research. 

In  recent  years,  significant  progress  in  modeling  the  structural  evolution  of 
materials  during  processing  has  been  made.  These  advances  have  been 
primarily  in  application  of  continuum  methods  where  phenomenological 
constitutive  relations  are  used  to  describe  how  a  material  will  behave. 
These  methods  can  be  used  to  predict  the  final  shape  of  a  specimen 
following  deformation  processing  and  its  temperature  history.  However, 
these  models  have  had  little  impact  on  material  design  because  they 
cannot  quantitatively  determine,  predict,  or  manipulate  the  internal 
structure  of  the  material.  Novel  materials  are  often  manufactured  via 
advanced  processing  techniques  that  greatly  influence  the  structure  of  the 
material  and  thus  its  mechanical  properties.  It  is  the  internal  structure 
that  controls  the  properties  of  a  material,  and  it  is  the  primary  mechanism 
that  materials  scientists  and  engineers  have  to  optimize  material 
properties.  Material  optimization  requires  researchers  to  develop  new 
material  models  that  link  theories  for  the  properties  and  performance  of 
materials  with  the  thermo-mechanical  processing  steps  required  to 
produce  them.  Ideally,  these  models  will  be  based  upon  understanding  the 
relationship  between  composition,  structure,  and  properties  and  the 
ability  to  process  materials  to  achieve  target  compositions  and  structure. 

Molecular  dynamics  simulations  are  used  to  determine  the  structural 
response  at  the  atomistic  scale  and  thus  provide  data  to  formulate 
constitutive  equations  for  the  coarse-grained  simulations.  For  example, 
the  structural  response  across  grain  boundaries  can  be  modeled  to  atomic 
fidelity  using  MD;  in  coarse-grained  simulations  the  boundary  would  be 
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modeled  as  simple  finite  elements  having  properties  derived  from  the 
small-scale  MD  simulation. 

MD  can  be  applied  to  understanding  synthesis  processes  to  some  extent, 
but  is  greatly  limited  in  this  application  because  of  the  time  scales 
involved.  To  meet  requirements  for  stability  and  accuracy,  the  time-step 
for  MD  is  small  relative  to  the  time  required  to  capture  diffusion  processes 
associated  with  synthesis  processes.  Also,  for  sintering  processes,  the 
length  scale  of  individual  grains  makes  the  application  of  MD  nearly 
intractable,  giving  rise  to  alternative  techniques  described  in  Chapter  5. 

Critical  to  accurate  MD  simulations  are  the  potentials  that  define  the 
forces  associated  with  atomic  interactions.  These  potentials  are  not 
fundamental  in  the  sense  that  the  potential  is  not  a  direct  result  of 
quantum  mechanical  computations.  Rather,  the  potentials  are  empirical 
relationships  that  are  inferred  from  a  combination  of  experimental  and 
computational  results.  It  is  common  to  base  the  potential  on  binary 
interactions  between  atoms,  although  it  is  often  necessary  to  consider 
multiple  bond  potentials.  The  selection  of  the  appropriate  potential  is 
central  to  the  science  of  MD;  no  general  atomic  potentials  fit  all  cases. 

A  review  of  literature  has  identified  several  technical  challenges  involved 
with  the  design  and  development  of  advanced  materials  at  the  nanoscale 
to  mesoscale: 

•  Develop  accurate  reactive  empirical  potentials. 

•  Develop  high-fidelity  3-D  digital  microstructure  models. 

•  Develop  methods  for  validating  the  digital  models. 

•  Develop  high-fidelity  atomistic  polycrystalline  models. 

•  Define  the  interface  between  the  atomistic  simulations  and  mesoscale 
simulations. 

•  Develop  methods  for  validating  atomistic  and  mesoscale  simulation 
results  with  experiment. 

2.1  Accurate  reactive  empirical  potentials 

The  computation  resources  required  to  run  ab  initio  MD  simulations  on 
systems  containing  millions  of  atoms  are  prohibitive.  Classical  potentials 
do  not  suffer  from  the  same  computational  restrictions  as  the  ab  initio  MD 
simulations.  However,  because  the  classical  potentials  are  only  model 
representations  of  interatomic  interactions,  classical  simulations  are 
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limited  by  the  accuracy  of  these  models.  For  these  models  to  be  useful, 
they  must  reproduce  as  closely  as  possible  the  experimental  or  ab  initio 
information  that  is  used  to  produce  them.  Attempts  have  been  made  to 
develop  empirical  potentials  that  are  “transferable,”  i.e.,  force  field 
parameters  fit  for  a  given  set  of  elements,  and  then  expect  the  same  force 
field  parameters  to  give  valid  results  for  completely  different  sets  of 
elements.  Another  approach  is  to  develop  a  systematic  parameter 
optimization  scheme  that  can  generate  the  force  fields  for  a  wide  variety  of 
molecular  systems. 

Reactive  empirical  potentials,  potentials  capable  of  making  and  breaking 
chemical  bonds  during  a  simulation,  are  required  for  molecular  dynamics 
simulation  of  chemical  reactions  that  take  place  during  the  sintering 
process.  Various  reactive  empirical  potentials  have  been  evaluated  to 
determine  their  suitability  for  modeling  the  sintering  of  SiC.  James  Lill 
(2010)  evaluated  the  Reactive  Empirical  Bond  Order  potentials  of  Brenner 
(Brenner  1990);  the  Adaptive  Intermolecular  Reactive  Empirical  Bond 
Order  (AIREBO)  and  the  potentials  of  Stuart  et  al.  (2000);  as  well  as  the 
Bond  Order  (BO)  potentials  of  Professor  David  Pettifor  (Pettifor  1989). 
Potentials  developed  by  Yu,  et  al.  (2007)  (Shan_20io)  and  the  Charge 
Optimized  Many  Body  (COMB)  potential  that  allows  for  dynamical 
electron  transfer  and  bonding  were  also  examined,  in  addition  to  the 
ReaxFF  formalism  developed  by  Professor  William  Goddard  and  Associate 
Professor  Adri  van  Duin.  ReaxFF  (van  Duin  et  al.  2001;  Chenoweth  et  al. 
2008)  is  perhaps  the  most  sophisticated  reactive  force  field  yet  proposed. 

Currently,  the  Tersoff  potential  in  Large-scale  Atomic/Molecular 
Massively  Parallel  Simulator  (LAMMPS)  (Plimpton  1995)  is  being  used  to 
run  the  atomistic  simulations  on  SiC.  Based  on  Lill’s  Reactive  Potential 
Survey,  the  project  “Using  AIREBO  Potentials  to  Optimize  Properties  of 
SiC  Composites  Reinforced  with  Carbon  Nanotubes”  (Lill  et  al.  2011)  was 
initiated  to  develop  an  optimized  parameter  set  for  modeling  the  sintering 
of  SiC  using  AIREBO  potential.  This  work  will  develop  new  AIREBO  force 
fields  to  support  ERDC’s  “Super  Carbon  Nanotube  Tensile  and 
Compressive  Infrastructure  Materials”  project  designed  to  develop  “super” 
composite  materials  with  dramatically  improved  compressive  strengths, 
tensile  strengths,  and  stiffness,  as  well  as  strength-to-weight  and  stiffness- 
to-weight  ratios.  Deliverables  will  include  the  following: 

•  Incorporation  of  new  AIREBO  force  fields  into  LAMMPS 
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•  Other  enhancements  to  LAMMPS  and  Vienna  Ab  initio  Simulation 
Package  (VASP)  including  an  option  to  perform  constant-temperature, 
constant-stress  (NST)  dynamics 

•  The  development  of  a  training  set  to  be  used  with  a  Multiple  Objective 
Evolution  Strategies  (MOES)  to  evolve  optimized  AIREBO  parameter 
sets 

•  Optimized  AIREBO  force  fields  for  the  SiC  +  CNT  system 


2.2  Interface  between  the  atomistic  simulations  and  mesoscale 
simulations 

The  validation  between  the  atomistic  and  mesoscale  simulations  has  not 
yet  started  in  the  ERDC  research  program.  A  target  system  size  has  been 
selected  with  which  to  conduct  the  validation.  A  system  containing  3  by  3 
by  3  10-mm  grains  each  was  selected  as  the  system  size  that  would  be 
modeled  for  both  the  atomistic  and  mesoscale  simulations  (Figure  1). 

2.3  Validating  atomistic  and  mesoscale  simulation  results  with 
experiment 

The  atomistic  simulations  are  progressing  to  the  point  of  coupling  them 
with  the  mesoscale  models.  The  experimental  work  is  progressing 
independent  of  the  simulations  for  now.  The  preliminary  validation 
between  the  atomistic  and  mesoscale  simulations  with  experimental 
results  has  started  (Figure  2)  and  will  continue  to  be  updated  as  the 
software  for  constructing  the  3-D  digital  models  is  progresses  and  the 
sintering  capabilities  have  been  developed. 

2.4  Example:  Effects  of  carbon  nanotubes  and  graphene 

SiC  is  a  remarkably  hard  and  strong  ceramic  material  consisting  of 
tetrahedral  structures  of  carbon  and  silicon.  Graphene  is  a  one-atom-thick 
planar  sheet  of  carbon  atoms  arranged  on  a  hexagonal  crystal  lattice.  Sp2 
bonds  connect  carbon  atoms  in  graphene,  which  are  the  strongest  bonds 
found  in  nature.  As  of  2009,  graphene  and  carbon  nanotubes  appear  the 
strongest  material  ever  tested,  with  tensile  strengths  of  about  105  GPa 
(15.5-million  psi),  or  150  times  that  of  high-strength  steel. 
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Atomistic  representation  of  SiC  Grain 


DEM  representation  of  SiC  Grain 


Embedded  C-allotrope 


Figure  1.  Discrete  element  models  (DEM)  of  reinforced  SiC  generated  to  match  atomistic 
structures  simulated  in  MD— grain  for  grain,  linker  for  linker.  MD  and  DEM  simulations  will  be 
run  for  comparative  analyses.  Results  from  MD  simulations  will  be  used  to  parameterize  DEMs. 


Figure  2.  Subsequent  DEM  simulations  for  reinforced  SiC  systems  at  larger  length  scales, 
comparable  with  those  of  the  experiment.  Comparative  studies  may  include  tests  of  hardness 

and  tensile  strength. 
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Results  from  preliminary  studies  carried  out  at  the  Information 
Technology  Laboratory,  Vicksburg,  MS,  of  the  U.S.  Army  Engineer 
Research  and  Development  Center  suggest  it  might  be  feasible  to 
synthesize  graphene-embedded  SiC  composites  with  deposition  methods. 
Classical  MD  simulations  were  performed  to  study  the  re-depositing 
sublimated  SiC  on  pristine  and  porous  sheets  of  graphene  (Figure  3). 


. . . 

v  t.  

- i _ t _ ^ _ ■ _ I 

K.  Suenaga/AIST 

Nanochains 


A.  Van  Duin  et  at. 

Nanorings  Nanoribbons 


Figure  3.  In  the  initial  setup,  a  pristine  graphene  sheet  comprising  480  graphene  atoms  was 
positioned  near  a  hot  gas  composed  of  2000  Si  and  C  atoms  in  a  240A  x  110A  x  100A  MD 
cell  with  periodic  boundary  conditions  on  all  sides.  In  the  second  simulation,  30  percent  of 
the  atoms  in  the  graphene  sheet  were  removed  to  create  vacancy  defects,  but  all  other  setup 

parameters  were  the  same. 


Each  simulation  time-step  (ts)  is  equal  to  1  femto  second  (fs)  in  real  time 
(1  ts  =  1  fs).  One  million  time-steps  are  required  to  simulate  1.0  ns  in  real 
time.  Each  simulation  ran  for  a  total  of  -1.2  ns.  Temperature  was 
maintained  between  1000  and  2000  K,  well  within  the  solid-phase  regime 
for  SiC.  Monoatomic  Si  and  C  were  observed  to  nucleate  at  the  edges  of 
pristine  graphene,  and  at  the  edges  and  vacancy  defects  of  porous 
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graphene.  Smaller  nanoparticles  rapidly  precipitated  (in  <  l  ns)  from  the 
gas  to  aggregate  around  the  graphene  as  well. 

It  is  expected  that  this  process  will  continue  and  that  solid  polycrystalline 
matrices  of  SiC  will  eventually  form  around  graphene,  wherein  strips  of 
graphene  will  interlace  through  multiple  crystalline  domains. 

It  is  hypothesized  that  when  compressed  and  sintered,  interlacing  and 
possibly  interconnecting  strips  of  graphene  will  reinforce  the  SiC  matrix  to 
increase  the  toughness  or  tensile  strength  of  the  overall  structure. 

2.5  Challenges 

Challenges  would  involve  optimizing  graphene  strip  width,  length, 
distribution,  and  interconnectivity  for  greater  tensile  strength  while 
avoiding  fragmentation  of  the  SiC  matrix;  and  optimizing  SiC  grain  size 
and  morphology  to  maintain  compressive  strength  and  hardness. 

Classical  MD  can  be  used  to  investigate  the  mechanical  properties  in 
reinforced  SiC  composites  at  the  nanoscale,  and  particularly  how  SiC 
might  be  engineered  at  the  nanoscale  to  resist  fracture. 

What  is  not  clear  is  how  the  presence  of  the  graphene  between  the  grains 
might  affect  the  structure  of  the  grain  boundary  region.  Typically,  the 
atomic  structure  within  the  grain  boundary  is  amorphous,  possibly 
because  it  represents  a  linkage  between  the  mismatch  in  crystalline 
structures  of  the  two  grains  making  up  the  boundary.  How  will  the 
presence  of  the  crystalline  graphene  within  the  normally  amorphous 
region  affect  this  structure?  Will  the  graphene  affect  the  strength  and 
toughness  of  the  grain  boundary  beyond  acting  as  a  simple  tensile 
element?  Atomic-scale  studies  can  complement  experimental  efforts  to 
answer  these  questions  (Walker  et  al.  2011). 

2.6  Progress 

Molecular  dynamics  simulations  have  been  performed  to  study  structural 
and  mechanical  properties  of  graphene-based  SiC  ceramics.  SiC  exhibits 
exceptional  compressive  strength  and  hardness,  but  poor  tensile  strength 
and  fracture  toughness.  Graphene  is  currently  accepted  as  the  strongest 
material  ever  tested,  with  tensile  strength  200  times  greater  than  steel. 
This  research  seeks  to  address  challenges  in  the  development  of 
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graphene-embedded  SiC  ceramic  composites  with  significantly  increased 
tensile  strength  and  fracture  toughness. 

Graphene  dissociates  at  2800  °C,  while  SiC  vaporizes  at  2750  °C.  It  is 
believed  that  within  this  temperature  range  particular  mixtures  of  gaseous 
SiC  and  graphene  can  combine,  condense,  and  be  annealed  into  graphene- 
interlaced  SiC  composites  with  enhanced  mechanical  properties.  The 
relationships  between  atomic  structure  and  mechanical  properties  need  to 
be  understood  to  determine  whether  interlacing  SiC  grains  with  C- 
allotropes  can  increase  tensile  strength,  similar  to  rebar  in  concrete  or 
straw  in  adobe.  How  hardness  and  compressive  strength  of  SiC  are 
affected  when  various  C-allotropes  are  embedded  needs  to  be  determined. 

Images  in  Figures  4  and  5  are  snapshots  from  some  of  the  first  simulations 
carried  out  to  investigate  the  mechanical  response  of  CNT-reinforced  SiC. 
They  show  two  slab-shaped  wurzite  crystal  grain  segments,  separated 
vertically  by  an  amorphous  layer  of  SiC,  interlinked  by  a  (26,0)  CNT  20A 
in  diameter  and  8A  in  length.  The  0001  directions  of  the  crystals  in  each  of 
the  grain  segments  are  aligned  for  this  case.  The  system  consists  of 
-50,000  Si  and  C  atoms.  Since  the  effect  of  Hall-Petch  hardening  is 
greatest  in  SiC  when  the  average  grain  diameter  is  between  5  nm  and 
20  nm,  the  dimensions  of  the  segments  were  chosen  to  explore  this  range. 


Figure  4.  Simulations  of  reinforced  SiC  grain  segment  failing  under  tension.  Snapshot 
sequences  show  MD  simulations  of  mechanical  response  of  CNT-reinforced  SiC  grain 
segments  to  tensile  stress.  The  purpose  is  to  determine  if  embedded  CNTs  can  increase 

tensile  strength. 
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Figure  5.  Simulations  of  reinforced  SiC  grain  segment  failing  under  compression.  Snapshot 
sequences  show  MD  simulations  of  mechanical  response  of  CNT-reinforced  SiC  grain  segments  to 
compressive  stress.  The  aim  is  to  observe  the  effect  of  the  embedded  CNT  on  hardness. 

Figure  4  shows  the  SiC-reinforced  segment  in  initial,  intermediate,  and 
final  points  in  time.  Separation  initiates  at  the  amorphous  buffer  and 
continues  across  the  amorphous  interface  until  the  two  segments  are 
separated.  Likewise,  Figure  5  shows  the  SiC-reinforced  segment  under 
compression  to  understand  how  the  inserted  CNT  affects  the  overall 
hardness  of  the  material.  Simulations  are  still  ongoing  to  study  the  effects 
of  CNTs  of  different  sizes  and  numbers  on  fracture  in  SiC. 

2.7  Technical  gaps 

Nanotechnology  is  a  relatively  new  discipline  in  materials  science.  Most  of 
the  technical  gaps  are  associated  with  the  ability  to  study  and  control 
materials  with  atomic  and  molecular  precision,  described  as  follows: 

•  Experimental  gaps  in  the  ability  to  characterize  the  atomic  structure 
and  composition  of  materials. 

•  Simulation  gaps  in  the  ability  to  construct  high-fidelity  3-D  digital 
models. 

•  Computational  gaps  in  the  reactive  empirical  potentials. 

•  Gaps  in  the  ability  to  conduct  verification  and  validation  on  all  levels. 

2.8  Future  plans 

The  work  described  in  this  chapter  is  still  in  progress,  and  it  will  take  most 
of  FY11  to  complete.  In  FY11,  simulations  using  the  models  will  begin  as 
they  become  available.  This  should  happen  toward  the  end  of  the  second 
quarter  of  FY11  with  the  simulation  of  the  interaction  between  particles  of 
differing  size,  grain-grain  interactions,  and  the  effect  of  grain  boundary 
size  and  composition  of  grain  growth  during  sintering.  Virtual  experiments 
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can  be  conducted  once  the  MBuilder  and  atomistic  model  results  are  avail¬ 
able.  Simulations  of  CNT  and  CNT  fibers  embedded  in  a  SiC  matrix  will  be 
the  last  thing  to  come  online,  and  simulations  of  these  systems  probably 
will  not  start  until  the  end  of  FY11. 

The  software  developed  during  FY11  will  form  the  basis  for  all  future  simu¬ 
lations.  This  software  must  be  robust  enough  to  capture  all  modeling  re¬ 
quirements  for  future  simulations.  The  objective  of  the  modeling  and  sim¬ 
ulations  is  to  support  experimental  efforts  to  understand  and  control  the 
sintering  process.  One  of  the  most  challenging  components  of  nanotech¬ 
nology  is  the  precise  control  of  the  atomic  and  molecular  structure  of  the 
material  during  processing.  The  same  level  of  precision  is  required  for 
characterizing  the  materials  and  construction  for  the  3-D  digital  represen¬ 
tation  of  the  materials  used  in  the  simulations. 
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3  Modeling  of  Polycrystalline  Materials 

The  mesoscale  represents  the  scale  that  is  larger  than  atomic  scale,  yet  like 
the  atomic  scale,  is  composed  of  interacting  discrete  components,  each 
consisting  of  millions  of  atoms.  In  many  respects,  mesoscale  mechanics  is 
similar  to  granular  mechanics  in  that  the  mechanical  response  is  as  much 
controlled  by  interactions  between  discrete  objects  as  by  the  atomic 
makeup  of  the  objects  themselves.  The  atomic  interaction  controls  how  the 
grains  are  formed  and  the  nature  of  the  grain  boundary;  but  once  formed, 
the  mechanical  response  depends  on  the  interaction  of  the  grains. 

To  design  materials  to  achieve  properties  for  specific  applications,  a  large 
number  of  parameters  have  to  be  evaluated.  An  experimental  approach 
based  on  trial-error  becomes  prohibitively  time-consuming  and  expensive, 
making  attractive  an  approach  based  on  numerical  simulations.  Numeri¬ 
cal  simulations  have  the  advantage  of  being  simpler  and  cheaper,  but 
provide  data  about  details  not  assessable  through  experimental  measure¬ 
ment.  Ideally,  MD  models  that  include  all  the  atoms  should  be  used  to 
capture  the  material  responses  fully.  However,  such  comprehensive  MD 
models  are  not  feasible  considering  the  expense  and  limitation  of  compu¬ 
tational  power.  As  an  alternative,  mesoscale  models  focus  on  the  size, 
shape,  orientation  and  spatial  arrangement  of  individual  components  and 
the  properties  and  characteristics  of  interaction  between  components  in 
lieu  of  modeling  all  atoms  explicitly. 

The  granular  nature  of  the  mesoscale  is  derived  from  the  scale  of  the  het¬ 
erogeneity  whereby  atomic  attractions  acting  at  the  boundaries  between 
discrete  components  are  much  smaller  than  those  between  atoms  within 
the  bulk  of  the  components.  In  a  general  sense,  therefore,  materials  at  the 
mesoscale  may  be  considered  to  be  large  assemblies  of  components  or 
grains  of  different  sizes,  shapes,  orientations,  and  compositions  that  are 
bonded  together  to  form  continuous  bulk  materials.  The  sizes  of  the  com¬ 
ponents  can  range  from  nanometers  to  millimeters.  The  properties  and 
performances  of  the  materials  result  from  an  amalgam  of  properties  and 
responses  of  individual  components  and,  equally  important,  interactions 
between  components.  The  size,  shape,  compositions,  orientation,  and 
spatial  arrangement  of  individual  components  and  the  properties  and 
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characteristics  of  interaction  between  components  determine  the  strength, 
fracture  toughness,  ductility,  and  fatigue  life  of  a  material  (Groeber  et  al. 
2008),  often  to  a  greater  extent  than  the  properties  of  the  components 
themselves.  In  that  sense,  mesoscale  has  a  great  deal  in  common  with 
granular  mechanics,  and  many  of  the  methods  used  to  model  granular 
media  are  equally  applicable  to  modeling  the  mesoscale. 

For  polycrystalline  materials,  the  Hall-Petch  effect,  in  which  the  yield 
strength  of  many  materials  increases  with  the  decrease  of  grain  sizes  by 
increasing  dislocation  pileups  at  the  grain  boundaries,  is  reversed  when 
the  sizes  of  grains  fall  into  the  10-  to  100-nm  range.  At  this  point  the 
deformation  mechanisms  change  from  dislocation  motions  across  grains 
to  grain  boundary  diffusion  (Schiotz  and  Jacobsen  2003).  A  notable 
example  is  the  superior  properties  of  biologic  materials  that  are  possible 
because  the  otherwise  weak  and  brittle  components  making  up  the 
materials  have  such  small  dimensions. 

Natural  materials  derive  their  superior  performance  by  incorporating 
small-scale  particles  into  structures  that  inherently  promote  toughness. 
Lee  et  al.  (2001)  show  that  elongated  grains  enhance  the  fracture 
toughness  of  SiC  ceramic  by  crack  bridging  and  crack  deflection.  Yang 
(2007)  shows  that  deformation  mechanisms  of  Face-Centered-Cubic  and 
Hexagonal-Close-Packing  metals  strongly  depend  on  the  grain  orienta¬ 
tions.  Natural  structural  materials  such  as  nacre,  bone,  fish  scale,  wood, 
bamboo,  deer  antler,  and  many  others  (Espinoza  et  al.  2009;  Bruet  et  al. 
2008;  Fratzl  2007;  Yao  et  al.  2010;  Rhee  et  al.  2009;  and  Trim  et  al.  2010) 
illustrate  the  importance  of  spatial  arrangements  of  individual  compo¬ 
nents  or  grains  because  their  extraordinary  properties  largely  depend  on 
their  hierarchical  arrangements  of  individual  components.  Natural 
structural  materials  are  the  best  examples  of  how  the  size,  shape, 
compositions,  orientation,  and  spatial  arrangement  of  individual 
components  and  the  properties  and  characteristics  of  interaction  between 
components  all  work  together  to  achieve  desired  mechanical  properties 
and  responses  of  bulk  materials. 

SiC  ceramic  is  a  polycrystalline  material  with  the  size  of  grains  ranging 
from  hundreds  of  nanometers  to  microns  (Cao  et  al.  1996).  The  crystalline 
structures  of  SiC  include  three  major  types:  6H(a),  4H(a),  and  3CO).  The 
grain  boundaries  of  SiC  have  thicknesses  of  a  few  nanometers  and  can  be 
amorphous  in  structure  largely  depending  on  the  sintering  conditions 
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(Tsurekawa  et  al.  1995).  As  mentioned  earlier,  SiC  has  high  compressive 
strength  but  limited  tensile  strength  and  fracture  toughness  because 
atomic  bonding  is  covalent.  The  failure  mechanism  of  SiC  bulk  materials 
can  be  cleavage  through  grains  (transgranular  failure)  or  cracks  along  the 
grain  boundaries  (intergranular  failure).  A  lot  of  research  has  been 
performed  to  increase  the  fracture  toughness  or  tensile  strength  of  the  SiC 
materials  (Shinozaki  et  al.  1993).  Cao  et  al.  (1996)  used  sintering  additives 
to  achieve  different  fracture  toughness  by  developing  microstructures  with 
different  grain  sizes  and  shapes.  Kim  et  al.  (1997)  also  used  large  ellipsoid¬ 
shaped  grains  among  smaller  and  more  equiaxial-shaped  grains  to 
improve  the  fracture  toughness  of  bulk  materials.  LaSalvia  et  al.  (2010) 
also  demonstrated  that  the  fracture  toughness  of  the  SiC  bulk  material  was 
improved  by  weakening  the  grain  boundaries  and  allowing  cracks  to  grow 
intergranularly  rather  than  transgranularly.  Zhou  et  al.  (2001)  enhanced 
Young’s  modulus,  hardness,  flexure  strength,  and  fracture  toughness  by 
changing  the  chemistry  of  intergranular  phase  in  the  SiC  ceramic  material. 
Fibers  have  been  incorporated  into  SiC  material  to  improve  the  bulk 
fracture  toughness  (Wang  et  al.  2006).  Carbon  nanotubes  with  their 
extraordinary  strength  and  properties  are  also  being  used  as 
reinforcement  in  SiC  (Karandikar  et  al.  2008).  Determination  of  howto 
take  advantage  of  the  properties  of  CNTs  to  improve  the  tensile  strength 
and  fracture  toughness  of  SiC  requires  detailed  modeling  at  multiple 
scales.  Numerical  simulation  of  SiC  materials  at  the  mesoscale  can  be  a 
powerful  tool  for  learning  how  to  configure  these  materials  for  high- 
strength  and  high-fracture  toughness. 

To  build  detailed  mesoscale  models,  the  detailed  mesostructure  of 
materials  needs  to  be  digitally  characterized  and  implemented  in  the 
computer  models.  To  date,  Electron  Back-Scatter  Diffraction  (EBSD)  serial 
sectioning  is  the  most  accurate  way  of  characterizing  samples  and  building 
digital  3-D  mesostructure  for  computer  models  (Rollett  et  al.  2007).  The 
mesostructure  can  also  be  built  synthetically  by  building  statistical 
representation  of  the  mesostructure  using  data  from  experiments  (Rollett 
et  al.  2007).  Model  construction  methods  are  critical  to  the  success  of 
mesoscale  modeling  and  are  described  in  detail  in  the  section  “High- 
fidelity  atomistic  polycrystalline  models.” 

After  mesostructured  detail  is  implemented  in  the  computer  models, 
suitable  modeling  techniques  are  needed  to  study  the  mechanical  prop¬ 
erties  and  deformation  mechanisms  of  the  mesostructure.  Researchers  are 


ERDC  TR-12-9 


22 


using  a  number  of  methods.  The  finite  element  (FE)  method  combined 
with  cohesive  elements  to  represent  fractures  is  considered  the  traditional 
approach  for  such  problems.  Meshless  approximations  such  as  the 
peridynamics  method  have  begun  to  garner  greater  interest  because  of 
their  perceived  ability  to  model  fracture  without  using  cohesive  elements 
with  their  predefined  orientations.  The  DEM  is  also  being  used  for 
problems  for  which  the  grains  have  sufficient  intergranular  mobility  to 
which  continuum  mechanics  do  not  strictly  apply. 

The  mesostructure  models  are  not  only  built  to  study  the  deformation 
mechanisms  within  the  mesoscale  but  also  to  generate  constitutive  laws  to 
be  used  for  macroscale  models.  The  question  thus  arises,  what  are  the 
smallest  mesoscale  models  that  can  adequately  capture  the  physics  in 
mesoscale  yet  define  a  sufficiently  large  domain  to  constitute  a  valid 
representative  elementary  volume?  Appropriate  homogenization 
techniques  have  to  be  examined. 

Structural  SiC  is  produced  by  initially  mixing  fine  (submicron)  and  pure 
SiC  powder  with  nonoxide  sintering  aids.  The  powdered  material  is 
compacted  and  is  then  sintered  under  high  temperature  and  sometimes 
under  high  pressure.  The  sintering  of  SiC  can  be  solid  phase  or  liquid 
phase  (Cao  et  al.  1996).  As  mentioned  previously,  to  achieve  desired 
mechanical  properties,  the  SiC  mesostructure  needs  to  be  controlled. 
Modeling  the  sintering  process  at  the  mesoscale  can  be  instrumental  in 
guiding  the  selections  of  powders,  sintering  temperature,  pressure,  and 
sintering  aid.  However,  the  physics  of  sintering  are  more  involved  than 
mechanical  response,  as  described  in  Chapter  5. 

In  the  following  sections  in  this  chapter,  different  mesoscale  modeling 
techniques  will  be  discussed  in  detail,  and  some  preliminary  simulations 
of  SiC  will  be  presented.  The  technical  gaps  in  this  area  will  also  be 
discussed. 

3.1  Mesoscale  modeling  methods 

3.1.1  Finite  Element  Method 

Espinosa  and  coworkers  (Zavattieri  and  Espinosa  2001,  2003;  Espinosa 
and  Lee  2002;  Espinosa  and  Zavattieri  2003)  performed  grain-level 
simulations  of  ceramic  materials  under  dynamic  loading.  The  grain 
structure  either  comes  from  digitalization  of  an  actual  microstructure  or  is 
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generated  from  Voronoi  tessellation.  Only  two-dimensional  grain 
structures  were  modeled.  The  grains  were  discretized  into  triangular  finite 
elements,  and  the  grain  boundaries  were  modeled  with  cohesive  elements 
inserted  between  grains.  The  elastic  and  thermal  anisotropic  model  is  used 
to  describe  the  single  crystal  anisotropic  behavior  stress  tensor  Sij  of  the 
grains: 


Sij  =  CijkiHki  (1) 

where  Hki  is  a  logarithmic  strain,  and  Qju  is  the  elastic  anisotropic 
material  stiffness  tensor  in  the  global  coordinates.  Each  grain  is  assumed 
to  be  elastic-orthotropic,  and  the  orientation  of  the  principal  material 
directions  differs  from  grain  to  grain.  The  cohesive  law  used  for  grain 
boundaries  is: 
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where 

X  =  nondimensional  effective  displacement  jump 
un,  ut  =  normal  and  tangential  displacement  jumps  at  the  interface 
8 n,  8f  =  critical  values  at  which  interface  failure  takes  place 

8 n  and  8f  depend  on  the  maximum  normal  traction  and  the  fracture 
toughness  at  the  interfaces. 

The  models  predefine  the  cracks  to  grow  along  grain  boundaries.  For 
dynamic  loading,  experimental  results  show  that  cracks  nucleate  at  triple 
joints  and  grow  along  grain  boundaries.  The  assumption  of  intergranular 
failure  of  models  is  justified  for  dynamic  loading.  However,  for  static 
loading,  Cao  et  al.  (1996)  and  Kim,  et  al.  (1997)  show  that  cracks  can  grow 
transgranularly  through  cleavage  for  brittle  ceramic  materials.  Current 
finite  element  models  cannot  capture  transgranular  crack  growth. 

Postek  et  al.  (2005)  also  used  the  finite  element  method  to  model  ceramic 
brittle  materials.  Their  models  considered  the  grain  boundaries  to  be  a 
thick  layer  of  materials  with  different  properties  from  the  materials  inside 
the  grains.  They  considered  the  material  inside  grains  to  be  elastic  and  the 
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grain  boundaries  to  be  elastic-plastic  materials  undergoing  finite  strain. 

The  failure  criteria  for  the  grain  boundary  material  were  based  on  the 
equivalent  plastic  strain.  Zhang,  et  al.  (2005)  also  used  the  finite  element 
method  to  study  3-D  Voronoi  polycrystal  (i.e.,  the  boundaries  of  the 
crystals  are  computer-generated  Voronoi  cells).  In  their  model,  crystal 
plasticity  laws  are  applied  to  the  grains,  and  cohesive  elements  are  used  to 
model  grain  boundaries. 

3.1.2  Peridynamics  method 

Peridynamic  theory  (Silling  and  Askari  2005;  Askari  et  al.  2008)  is 
employed  as  a  nonlocal  meshless  method  that  adopts  a  computational 
framework  similar  to  MD.  Unlike  classic  continuum  mechanics,  the  local 
balance  of  linear  momentum  is  expressed  as: 

pil(x,t)  =  f  f\u(X',t)-u(x,t),X' -X]dVx,  +b(X,t)  (3) 

hx 

where 

p  =  mass  density  in  the  reference  configuration 
u  =  displacement  vector  field 
x  =  position  of  particle  X 
x'  =  position  of  particle  X’ 

t  =  time 

Hx  =  a  neighborhood  volume  of  particle  X 
dVx ,  =  infinite  small  volume  around  particle  X’ 

f  =  a  pairwise  force  function  whose  value  is  the  force  vector  (per  unit 
volume  squared)  that  the  particle  X'  exerts  on  the  particle  X) 
b  =  body  force  density  field 

The  pairwise  force  function  includes  the  bonding  laws  between  material 
points.  The  peridynamic  method  is  implemented  in  a  general  molecular 
dynamics  code  LAMMPS  (Parks  et  al.  2010).  In  LAMMPS,  only  a 
microelastic  brittle  bonding  law  is  currently  implemented.  The 
microelastic  bonding  law  is  as  follows: 

f[u(x',0  -  u(x,t),x'  -  x]  =  C  J[u(x',f)  -  u(x,0]  -  (x'  -  x)||  (4) 

x'-x  "L  "  u(x',f)-u(x,f) 
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where  C  is  the  spring  constant  of  the  bond.  The  stretch  of  the  bond  is 
defined  as 


u(x',t) — u(x,t)l  +  (x/ — x) 


When  the  stretch  of  a  bond  reaches  the  critical  value,  the  bond  is 
considered  to  be  broken.  Damage  is  calculated  by  dividing  the  number  of 
broken  bonds  by  the  total  number  of  bonds. 

Silling  and  Askari  (2005)  used  the  peridynamic  method  to  study 
polycrystalline  fracture  of  silicon  in  two  dimensions.  They  generated 
Voronoi  polycrystals  with  a  grain-size  distribution  close  to  a  Weibull 
distribution,  with  each  grain  assigned  a  random  orientation.  The 
microstiffness  of  the  grains  depends  on  the  orientation  of  the  grains  and 
the  direction  of  the  bonds.  The  critical  bond  S0  at  which  a  bond  breaks  is 
associated  with  the  fracture  energy  and  the  peridynamic  horizon,  which  is 
the  cutoff  distance  of  the  bonds. 


5 


o 


(6) 


where 

Go  =  fracture  energy  of  the  material 
k  =  bulk  modulus 

8  =  horizon  (the  cut-off  radius  of  material  volume  around  particle 

A) 

The  critical  bonds  between  grains  are  related  to  grain  boundaries  as  fol¬ 
lows: 


where  p  is  the  interface  strength  coefficient. 

The  authors  attempted  to  use  the  peridynamic  method  to  study  a  6-grain 
SiC  polycrystalline  model.  The  LAMMPS  formulation  was  used  with  a 
simple  microelastic  brittle  bonding  law.  Initially,  a  3-D  6-grain 
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microstructure  with  a  bimodal  distribution  was  generated  using  Mcbuilder 
developed  by  Rollett  et  al.  (2007).  Based  on  Silling  and  Askari  (2005)  the 
following  equation  was  used  to  calculate  the  spring  constant  of  the  grains: 

C~  (8) 

no 

Equations  6  and  7  were  used  to  calculate  the  critical  bond.  The  bulk  elastic 
modulus  and  fracture  energy  are  from  Li  and  Bradt  (1987)  and  Kikuchi  et 
al.  (2005).  The  elastic  stiffness  of  the  bonds  between  the  grains  was 
determined  using  Equation  8.  A  preliminary  study  was  performed  on  the 
effects  of  relative  elastic  stiffness  and  critical  bond  between  inside  grains 
and  across  the  grains  (grain  boundaries).  \A  uniaxial  strain  was  applied 
on  the  specimen. 

Figure  6  shows  the  plot  of  a  peridynamic  model  of  a  six-grain  polycrystal. 
Figure  7  shows  the  plots  of  the  fracture  path  and  damage  when  the  stiff¬ 
ness  between  grains  is  half  of  that  inside  grains,  and  the  critical  stretch  be¬ 
tween  grains  is  the  same  as  that  inside  grains.  The  simulation  results  in 
Figure  8  show  that  cracks  grow  across  grains,  when  the  stiffness  of  the 
bonds  between  grains  is  twice  as  that  inside  the  grains  and  the  critical 
stretch  between  grains  is  the  same  as  that  inside  the  grains.  The  simula¬ 
tion  results  in  Figure  9  show  that  cracks  grow  through  the  middle  of  the 
specimen  as  primarily  transgranular  fractures  when  the  stiffness  between 
grains  is  the  same  as  that  inside  grains,  and  the  critical  stretch  of  the  bond 
between  grains  is  half  of  that  inside  grains.  The  simulation  results  in  Fig¬ 
ure  10  show  that  cracks  develop  at  random  locations,  when  the  stiffness 
between  the  grains  is  the  same  as  that  inside  the  grains,  and  the  critical 
stretch  between  grains  is  twice  that  inside  the  grains.  Note  that  the  grain 
orientations  are  not  considered  and  the  average  size  of  the  grains  is  in  mil¬ 
limeters,  which  is  not  realistic  for  SiC  grains.  However,  those  models  show 
that  both  elastic  properties  and  fracture  properties  of  grains  and  grain 
boundaries  affect  the  fracture  mechanisms  of  the  grain  structures. 

3.1.3  Discrete  Element  Method 

The  discrete  element  method  treats  a  material  as  an  assembly  of  inter¬ 
acting  particles  that  are  rigid  except  at  interparticle  contacts.  Steinhauser 
et  al.  (2008)  used  the  discrete  element  method  to  model  alumina  and  SiC 
under  loading  at  high  strain  rate  at  coarse-grain  level.  The  particles  are 
arranged  in  hexagonal  packing  within  each  grain,  whereas  the  particles  at 
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Figure  6.  Peridynamic  model  of  a  six-grain  polycrystal. 
The  colors  are  applied  to  distinguish  the  grains. 


(a)  Fracture  path.  The  colors  represent  different  grains. 


(b)  Damage.  The  colors  denote  the  extent  of  damage.  The 
red  denotes  that  damage  reaches  one.  When  damage 
reaches  one,  the  crack  develops. 


Figure  7.  Peridynamic  simulation  results.  The  stiffness 
between  grains  is  half  that  inside  the  grains.  The  critical 
stretch  between  the  grains  is  the  same  as  that  inside  the 
grains. 
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a.  Fracture  path.  The  colors  represent  different  grains. 
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b.  Damage.  The  colors  denote  the  extent  of  damage.  The 
red  denotes  that  damage  reaches  one.  When  damage 
reaches  one,  the  crack  develops. 

Figure  8.  Peridynamic  simulation  results.  The  stiffness 
between  grains  is  twice  that  inside  the  grains.  The  critical 
stretch  between  grains  is  the  same  as  that  inside  the  grains. 

grain  boundaries  are  packed  randomly.  The  bonding  forces  between 
particles  include  contact  forces  and  cohesive  bonding  force.  This  method  is 
able  to  capture  fracture  at  the  grain  boundaries  and  cleavage  across  grains. 

To  study  whether  DEM  methods  can  capture  possible  plastic  deformation 
in  the  grains,  small-scale  DEM  simulations  were  performed  to  study  how 
the  packing  arrangement  of  the  particles  may  affect  the  fracture  mecha¬ 
nisms.  In  this  model,  the  specimen  size  is  4onm  x  4onm  x  160  nm.  The 
applied  interparticle  forces  include  contact  normal  force,  friction  force, 
and  cohesive  normal  bonding  force.  The  bottom  of  the  specimen  is  fixed, 
and  a  downward  velocity  is  applied  on  the  top.  Figure  li  shows  the  defor¬ 
mation  and  internal  stress  in  the  vertical  direction  of  a  specimen  with 
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a.  Fracture  path.  The  colors  represent  different  grains. 
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b.  Damage.  The  colors  denote  the  extent  of  damage.  The 
red  denotes  that  damage  reaches  one.  When  damage 
reaches  one,  the  crack  develops. 

Figure  9.  Peridynamic  simulation  results.  The  stiffness 
between  grains  is  the  same  as  that  inside  the  grains.  The 
critical  stretch  between  grains  is  one-half  that  inside  the  grains. 

close  random  packing.  The  color  denotes  the  internal  stress  in  the  vertical 
direction.  Figure  n  shows  that  the  particles  form  force  chains— signified 
by  linear  groups  of  red  particles— under  compressive  loading.  Figure  12 
shows  the  deformation  and  internal  stress  in  the  vertical  direction  of  a 
specimen  with  a  face-centered  close  dense  packing.  Figure  12  shows  that 
the  specimen  slips  at  approximately  a  45-degree  angle  under  compressive 
loading.  However,  if  the  cohesive  force  is  high  (twice  that  in  Figure  12), 
Figure  13  shows  that  there  are  slip  bands,  but  the  specimen  fails  by  split¬ 
ting  in  the  middle.  Figures  11-13  show  that  packing  particle  in  Face- 
Centered  Cubic  (FCC)  packing  may  capture  plastic  deformation  such  as 
slip  bands,  but  the  interparticle  force  needs  to  be  tuned  to  capture  experi¬ 
mentally  observed  deformation  mechanisms  inside  the  grains.  Note  that 
the  DEM  particles  are  being  used  here  as  surrogates  to  finite  elements  as  a 
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means  to  accommodate  fracture  paths  which  are  not  preselected  as  in  the 
case  of  using  cohesive  elements.  Arguably,  such  a  model  lacks  many  physi¬ 
cal  aspects  of  the  real  material,  in  particular  an  accurate  depiction  of  wave 
propagation  within  the  grains.  These  simulations  do  demonstrate  that 
stiffness  as  well  as  strength  controls  qualitative  features  of  failure. 


a.  Fracture  path.  The  colors  represent  different  grains. 


b.  Damage.  The  colors  denote  the  extent  of  damage.  The 
red  denotes  that  damage  reaches  one.  When  damage 
reaches  one,  the  crack  develops. 


Figure  10.  Peridynamic  simulation  results.  The  stiffness 
between  the  grains  is  the  same  as  that  inside  the  grains.  The 
critical  stretch  between  grains  is  twice  that  inside  grains. 
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Figure  11.  Deformation  and  internal  stress  in  the  vertical 
direction  of  a  specimen  with  close  random  packing.  The 
color  denotes  the  continuum  stress  at  the  vertical  direction. 
Red  denotes  the  highest  stress. 


ParbclesnpLtnnge 
-1  7108-001 

1  282e-00 1 
8  536e-002 
■4  257e-002 

2  270e-004 


8 


a.  The  slip  just  initiated. 
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b.  The  slip  progressed  until  complete  failure.  The  color 
denotes  the  continuum  stress  at  the  vertical  direction.  Red 
denotes  the  highest  stress. 


Figure  12.  The  specimen  has  close,  dense  packing,  and 
cohesive  bonding  force  is  small. 
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Figure  13.  All  the  parameters  are  the  same  as  for  the 
specimen  shown  in  Figure  12  except  that  the  cohesive 
stress  is  twice  that  applied  in  the  specimen  shown  in 
Figure  12.  The  color  denotes  the  continuum  stress  at  the 
vertical  direction.  Red  denotes  the  highest  stress. 


3.2  Technical  gaps  in  mesoscale  modeling 

Based  on  the  information  provided  in  the  previous  sections,  the  following 
shortcomings  were  found  for  all  the  methods  described: 


•  The  finite  element  method  is  the  best  established  for  modeling  poly¬ 
crystalline  material  at  the  mesoscale.  This  method  can  address  possible 
plastic  deformation  inside  the  grains.  However,  it  can  address  the 
intergranular  cracking  only  where  the  orientation  of  the  potential 
fracture  is  known.  Although  cohesive  elements  can  be  placed  between 
boundaries  of  the  elements,  doing  so  makes  the  discretization,  which 
should  be  arbitrary,  a  property  that  controls  fracture  toughness.  Pre¬ 
diction  of  general  fractures  within  a  solid  medium  is  still  beyond  the 
capabilities  of  continuum  theories. 

•  Methods  based  on  peridynamics  can  address  failure  inside  the  grains 
and  at  the  grain  boundaries.  However,  it  does  not  address  possible 
plastic  deformations  inside  the  grains.  Moreover,  many  details  about 
the  actual  physical  mechanisms  captured  by  the  peridynamic 
formulation  have  yet  to  be  resolved. 

•  The  DEM  may  be  used  to  address  both  intergranular  and  transgranular 
cracking,  and  it  may  even  address  plastic  deformations.  As  noted 
previously,  the  method  is  not  fundamental  for  modeling  continuum 
media,  although  it  could  be  argued  that  the  method  is  no  less  suited  for 
this  application  than  peridynamics.  In  any  case,  for  the  method  to  give 
realistic  quantitative  results,  suitable  particle  packing  for  the  grains 
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and  grain  boundaries  and  appropriate  interparticle  forces  between 
particles  need  to  be  developed. 

3.3  High-fidelity  3-D  digital  microstructure  models 

Realistic  microstructure  and  nanostructure  models  for  finite  element, 
mesoscale,  or  atomistic  simulations  of  the  sintering  process  are  difficult  to 
produce  without  special-purpose  preprocessing  software.  Professor  An¬ 
thony  Rollett  at  Carnegie  Mellon  University,  Pittsburgh,  PA,  has  developed 
a  method  for  constructing  high-fidelity  3-D  digital  material  models.  ERDC 
initiated  a  project,  “High-fidelity  3-D  digital  polycrystalline  microstructure 
models,”  to  produce  the  models  needed  for  the  submicron-scale  simula¬ 
tions  of  the  sintering  process  of  SiC.  The  project  will  produce  models  that 
allow  the  evolution  of  material  structure  during  processing  to  be  studied  at 
resolutions  ranging  from  atomistic  to  microscale.  The  project  is  supported 
through  the  PETTT  program.  The  principal  investigators  on  the  project  are 
Anthony  Rollett  and  James  Larentzos,  High  Performance  Technologies, 
Inc.  (HPTi).  This  project  will  develop  the  technology  needed  to  produce 
accurate  3-D  digital  microstructure  models  needed  to  assist  experimental 
efforts  to  design  and  test  novel  high-strength  materials.  Polycrystalline 
materials  such  as  ceramics,  concrete,  and  some  biomaterials  include 
grains,  particles,  flaws,  and  reinforcements  with  different  morphologies. 
The  new  tools  will  be  used  to  build  accurate  digital  microstructures  with 
grains  with  large  aspect  ratios,  flaws,  and  fiber  reinforcement  and  will  be 
based  on  Carnegie  Mellon’s  Microstructure  Builder  (MBuilder)  tool.  The 
digital  microstructure  will  be  validated  against  the  microstructure  data 
collected  via  experimental  techniques.  The  digital  microstructures  will 
then  be  incorporated  into  atomistic,  peridynamics,  finite  element  numeri¬ 
cal  models  to  study  the  response  of  materials  at  the  multiple  scales.  The 
project  will  focus  on  the  generation  of  digital  macrostructures  with  flaws, 
grains  with  large  aspect  ratios,  irregularities,  and  fiber  reinforcements 
such  as  carbon  nanotubes  and  carbon  nanotube-based  fibers. 

Any  attempt  to  simulate  the  structural  evolution  of  materials  accurately 
during  processing  must  begin  with  a  faithful  digital  representation  of  the 
initial  state  of  the  material.  However,  capturing  the  initial  state  of  the 
complex  structural  geometry  of  a  material  remains  problematic.  For 
modeling  and  simulation  to  be  viable  tools  for  studying  the  behavior  and 
performance  of  materials,  research  scientists  must  provide  a  solid  link 
between  the  experimental  material  under  study  and  the  abstract  digital 
representation  of  the  material  used  in  the  simulations.  The  experimental 
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characterization  of  the  material  provides  this  link.  The  material  building 
process  employs  a  combination  of  experimental  and  computational 
techniques  to  capture  the  complex  structure  of  the  material.  This 
technique  then  creates  a  statistically  equivalent  representative  3-D  digital 
model  for  input  into  computer  simulations.  It  takes  geometric  and 
crystallographic  observations  from  two  orthogonal  sections  of  a 
polycrystalline  material  (Figure  14).  The  experimental  data  collected  from 
the  two  observation  planes  of  the  sample  material  are  used  to  construct 
experimental  distribution  functions.  The  distribution  functions  created 
during  the  scanning  process  are  used  to  guide  the  construction  of  the  3-D 
digital  structure  (Figure  15).  In  the  end,  the  characterization  criteria 
provide  a  quantitative  measure  for  gauging  the  accuracy  and  quality  of  the 
abstract  model.  Furthermore,  the  method  can  be  used  to  generate  a  wide 
variety  of  hypothetical  material  states  systematically  for  3-D  digital 
structures  for  computer  simulation.  These  hypothetical  models  provided  a 
practical  means  of  conducting  virtual  experiments  to  examine  material 
synthesis  under  a  wide  variety  of  experimental  conditions.  These  models 
will  have  the  predictive  capability  necessary  to  model  the  evolution  of  the 
material  structure  during  processing. 

3.4  Methods  for  validating  the  digital  models 

Three  critical  elements  of  this  technology  are  the  following: 

•  An  accurate  3-D  digital  representation  of  the  microstructure  of  the 
material 

•  An  efficient  method  for  generating  the  finite  element  mesh  needed  for 
modeling  and  simulations 

•  A  software  package  that  automates  the  process  as  much  as  possible. 


Figure  14.  The  material  is  scanned  using  an  automated  EBSD  in  a  Scanning  Electron 
Microscope  (SEM),  also  known  as  Orientation  Imaging  Microscopy  (OIM)  to  examine  two 

surfaces  of  the  material. 
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Figure  15.  Generation  of  geometric  grain  configuration  by  assigning  Voronoi  tessellations  to 

ellipsoids. 


The  first  step  in  the  process  is  to  collect  the  data  from  the  two  orthogonal 
observation  planes  of  a  polycrystalline  material.  The  MBuilder  uses  auto¬ 
mated  EBSD  in  a  Scanning  Electron  Microscope  (SEM),  also  known  as 
Orientation  Imaging  Microscopy  (OIM),  to  examine  two  surfaces  of  the 
material.  The  experimental  distribution  functions  are  obtained  from  the 
data  collected  from  the  orthogonal  section  planes  for  each  of  the  measured 
parameters.  At  present,  distribution  functions  are  generated  for  the  size, 
shape,  and  orientation  of  the  grains  and  the  crystallographic  orientations 
and  misorientations  across  grain  boundaries  for  each  grain  in  the  obser¬ 
vation  planes.  Here,  it  is  assumed  that  the  region  covered  by  the 
observation  planes  is  adequate  to  serve  as  a  representative  volume 
element  of  the  material  and  that  the  experimental  distribution  functions 
obtained  from  this  data  can  adequately  characterize  any  anisotropy 
associated  with  the  experimental  data.  To  simulate  the  microstructure 
geometry,  ellipsoids  (grains)  were  generated  according  to  the  experimental 
distribution  function.  Ellipsoids  were  chosen  to  represent  grain  geometries 
because  they  are  simple  geometric  forms  that  describe  both  the  size  and 
aspect  ratio  of  the  different  grains  that  compose  a  polycrystalline 
microstructure.  Grain  geometries  were  then  constructed  based  on  the 
packed  ellipsoids  using  a  continuous  nucleation  and  growth  model  to  fill 
space  (Rollett  et  al.  2004). 

Once  the  geometric  3-D  description  of  a  polycrystalline  microstructure  is 
complete,  the  crystallographic  information  is  assigned  to  each  grain 
forming  the  abstract  microstructure  such  that  the  distribution  of  grain 
orientations  is  equivalent  to  the  measured  distribution  in  the  real 
polycrystalline  material.  The  distribution  in  the  misorientations  across 
grain  boundaries  is  reproduced  using  a  Monte  Carlo  simulation  technique. 
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The  Monte  Carlo  simulation  uses  the  root-mean-squared  difference 
between  the  distribution  functions  in  the  real  and  simulation  materials  for 
the  acceptance  criteria  in  the  Monte  Carlo  simulation.  Changing  or 
swapping  the  orientations  of  the  grains  generates  trial  crystallographic 
orientations.  The  difference  in  the  distribution  functions  is  recalculated.  If 
the  trial  configuration  results  in  a  lower  difference  function,  the  trial 
configuration  is  accepted.  If  the  difference  function  is  greater  for  the  trial 
configuration  but  the  difference  is  small  enough,  it  is  probabilistically 
accepted.  Otherwise,  it  is  rejected.  This  process  continues  until  some 
convergence  criteria  (input  parameter)  are  met  at  which  point  the  process 
stops.  The  crystallographic  orientations  of  the  real  material  and  abstract 
model  are  considered  statistically  equivalent  to  the  extent  that  the 
differences  in  the  distribution  functions  have  been  minimized.  Once  the 
grains  are  produced,  a  finite  element  mesh  is  generated  for  each  of  the 
grains  (Figure  16)  (Cornwell  and  Noack,  2005). 


Figure  16.  Finite  element  mesh  generated  from  a  sample  of  aluminum  (the  model 
consists  of  134  randomly  colored  grains).  The  technique  can  be  used  to  produce 
microstructure  models  for  finite  element  simulations  and  nanostructure  models 
of  SiC  for  atomistic  simulations. 

3.5  High-fidelity  atomistic  polycrystalline  models 

The  MBuilder  project  will  capture  the  structure  of  the  materials  and 
provide  a  way  to  validate  the  structure  of  the  models  against  experimental 
data.  The  atomistic  simulations  will  require  a  way  to  populate  the 
structure  with  atoms  having  correct  crystal  structure  and  orientations. 
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This  information  is  also  provided  by  the  experimental  data.  Software  is 
being  developed  that  will: 

•  Generate  atomic  crystal  grains 

•  Control  crystallographic  orientation  for  each  grain 

•  Produce  correct  mismatch  distribution  at  grain  boundary 

•  Provide  a  method  for  controlling  the  size  and  composition  of  grain 
boundaries 

•  Construct  boundary  conditions  for  simulation 

•  Insert  CNTs  and  CNT  fibers  into  the  SiC  matrix. 

This  work  is  being  done  in  parallel  with  the  MBuilder  project.  The 
software  package  QHULL  is  used  to  produce  the  grain  structure  needed  to 
develop  the  atomistic  models  of  the  material  structures.  Initially,  simple 
grain  structures  are  used  to  develop  the  code  for  the  atomic  crystal 
structures.  This  technique  can  be  applied  to  any  closed  volume  and  will 
work  equally  well  with  the  data  that  are  provided  by  the  MBuilder  project. 
Figures  17  and  18  show  examples  of  systems  with  two  grains  and  a  grain 
boundary  with  different  orientations.  The  crystallographic  orientation  of 
each  grain  is  the  same,  but  the  software  will  be  extended  to  change  the 
crystallographic  orientation  of  the  grains  to  produce  a  crystallographic 
orientation  mismatch  between  adjacent  grains. 

Figure  19  is  a  system  with  nine  grains:  eight  grains  on  the  corners  and  one 
diamond-shaped  grain  located  at  the  center. 

These  simple  models  can  be  used  to  capture  the  effect  of  particle  size  on 
grain  growth  during  sintering.  They  can  also  model  the  interaction  at  the 
grain  boundaries  and  the  effect  of  grain  size  and  composition  on  grain 
growth.  Some  evidence  in  the  literature  suggests  that  MD  simulations  are 
able  to  capture  aspects  of  particle-particle  interactions  that  are  not  possi¬ 
ble  in  the  current  versions  of  mesoscale  and  continuum  models.  It  may  be 
possible  to  improve  the  coarser  grained  models  based  on  information  ob¬ 
tained  from  the  MD  simulations.  On  the  other  hand,  there  always  may  be 
information  about  the  sintering  process  that  is  obtainable  only  through 
atomistic  simulations. 

CNTs  and  CNT  fibers  have  not  begun  to  be  inserted  into  the  SiC  matrix 
yet.  The  software  for  constructing  the  CNTs  and  the  CNT  fibers  is  available 
from  a  previous  project  designed  to  study  the  design  of  CNT  fibers  and 
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Figure  17.  SiC  system  with  two  grains.  Grain  boundary  is 
1  nm  wide  and  aligned  perpendicular  to  the  z-axis.  The 
system  dimensions  are  10  x  10  x  10  nm.  Both  grains 
have  the  same  crystallographic  orientation. 


Figure  18.  SiC  system  with  two  grains.  Grain  boundary 
is  1  nm  wide  and  the  x-axis  is  rotated  15  deg  out  of  the 
xy-plane.  The  system  dimensions  are  10  x  10  x  10  nm. 
Both  grains  have  the  same  crystallographic  orientation. 


ERDC  TR-12-9 


39 


Figure  19.  SiC  system  with  nine  grains:  eight  grains  on 
the  corners  and  one  diamond-shaped  grain  located  at  the 
center.  Grain  boundaries  are  1  nm  wide.  The  system 
dimensions  are  10  x  10  x  10  nm.  All  of  the  grains  have 
the  same  crystallographic  orientation. 


membranes.  This  work  will  be  completed  in  FYli.  Also,  the  models  are 
constructed  using  a  single  crystal  structure  of  SiC.  The  modular  form  of 
the  software  will  allow  other  crystal  structures  to  be  added  on  an  as- 
needed  basis.  Once  the  MBulider  project  provides  the  structure  for  the  ma¬ 
terials,  simulations  can  begin  on  3-D  digital  material  models  constructed 
from  experimental  data.  At  this  point,  the  experimental  data  will  begin  to 
be  used  to  validate  the  simulation  results  with  the  results  from  experi¬ 
ments.  The  results  from  atomistic  and  mesoscale  simulations  will  be  able 
to  be  validated  with  experiment. 
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4  Material  Synthesis 

The  initial  nanoscale  studies  of  the  synthesis  of  ceramic  (polycrystalline) 
materials  involved  a  number  of  tasks.  A  critical  task  was  to  identify  the 
most  promising  sintering  methods.  Attendant  to  those  tasks  was  the 
acquisition  of  initial  ceramic  synthesis  equipment.  In  addition,  ceramic 
synthesis  experiments  were  required  to  establish  baseline  methods.  These 
experiments  involved  SiC  and  geopolymers  with  iterations  and 
refinements  of  the  simulation-synthesis  experimental  cycle. 

SiC  is  a  well-studied  class  of  materials  consisting  of  a  silicon  and  carbon 
crystal  structure  with  over  250  polytypes.  Since  SiC  is  not  abundant  in 
nature,  it  is  usually  grown  in  one  of  several  ways.  Single  crystals  of  low 
quality  can  be  made  using  the  Acheson  Process.  This  is  a  commercial 
process  that  involves  high  temperatures  and  long  processing  times.  Bulk 
single  crystals  of  high  quality  can  be  produced  using  the  Lely  Method  and 
the  Modified  Lely  Method.  Dense  polycrystalline  ceramics  are  formed  by 
compacting  and  sintering  small  SiC  grains  (powder  sintering). 

4.1  Polycrystalline  ceramics 

Much  processing  is  done  to  form  the  desired  starting  material  for 
sintering.  The  smaller  the  starting  powder,  the  more  likely  it  is  to  sinter. 
The  powder  is  compacted  and  then  heated  (sometimes  while  under 
pressure)  to  form  a  dense  body.  The  dense  body  comes  about  through  the 
diffusion  of  matter  between  adjacent  particles  to  form  what  are  called 
grains.  The  driving  principle  behind  the  sintering  process  is  the  tendency 
to  minimize  surface  energy  by  minimizing  surface  area  and  maximizing 
radius  of  curvature.  Between  the  grains  are  pores.  To  eliminate  the  pores, 
the  grains  are  allowed  to  grow.  However,  larger  grains  have  less  desirable 
mechanical  properties.  To  minimize  grain  growth,  pressure  can  be  used 
while  sintering  to  cause  the  grains  to  rearrange  themselves,  thereby 
eliminating  the  pores  without  enlarging  the  grains.  Another  method  used 
to  limit  grain  growth  is  to  increase  the  heating  rate  while  sintering.  This 
allows  the  diffusion  to  occur  before  grain  growth  takes  place.  Pressure  in 
addition  to  high  heating  rates  is  beneficial  for  reducing  void  and  particle 


sizes. 
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4.2  Sintering  methods 

Many  methods  of  sintering  SiC  have  been  identified,  although  the  common 
feature  of  all  is  an  SiC  powder  “compact”  that  is  heated  in  a  resistance 
furnace.  Densification  occurs  in  several  stages.  First,  the  particles  join  at 
their  contact  points  via  diffusion  driven  by  a  minimization  in  surface 
energy.  Continued  diffusion  produces  a  3-D  network  of  solid  matter  as 
well  as  a  3-D  network  of  open  space  (pores).  Grain  growth  causes  the 
porosity  to  decrease  until  the  pores  become  isolated  and  sometimes 
coalesce.  Densification  stops,  but  the  grains  may  continue  to  grow, 
sometimes  engulfing  the  pores. 

Several  factors  influence  the  ability  of  SiC  to  sinter.  Sintering  is  a 
competition  between  minimizing  surface  energy  of  the  particles  and  the 
grain  boundaries.  Smaller  particles  will  have  a  higher  incentive  to  sinter. 
However,  the  decrease  in  surface  energy  is  countered  by  an  increase  in 
energy  of  the  grain  boundaries.  Another  factor  that  affects  the  final 
sintered  density  is  the  green  density  (i.e.,  the  density  of  the  power  before 
final  sintering).  Usually,  the  higher  the  green  density,  the  higher  will  be 
the  sintered  density. 

Either  solid-state  methods  or  liquid-state  methods  can  achieve  sintering. 
Liquid-state  sintering  is  similar  to  solid-state  sintering  except  for  the 
addition  of  sintering  additives  to  the  SiC  powder.  Typically,  these  are 
metal  oxides  (yttria,  alumina,  magnesia).  These  additives  form  a  eutectic 
composition  and  melt  at  the  sintering  temperature  forming  a  liquid  phase 
between  the  SiC  grains.  This  liquid  effectively  reduces  the  friction  between 
the  SiC  particles,  allowing  them  to  move  around  and  settle,  and  increases 
the  rate  of  diffusion.  This  procedure  is  typically  done  in  conjunction  with 
pressureless  sintering,  hot  pressing,  or  sintering  forging.  The  only 
equipment  needed  is  that  for  solid-state  sintering.  The  advantages  of 
liquid-state  sintering  are  lower  sintering  temperatures,  higher  density,  and 
less  grain  growth.  Nevertheless,  since  the  liquid  phase  is  located  in  the 
grain  boundary,  the  strength  of  this  phase  is  a  limiting  factor  to  the 
strength  of  the  polycrystalline  composite.  In  addition,  the  working 
temperature  range  cannot  exceed  that  of  the  liquid  phase. 

In  hot  pressing,  heat  and  pressure  are  applied  to  a  sample  (up  to  2000  °C 
and  50  MPa)  and  can  be  done  in  a  vacuum.  Hot  pressing  is  less  effective 
for  fine  powder  samples  that  are  already  well  compacted  because  the 
compaction  process  has  already  reduced  the  void  space.  Special  hot  press 
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equipment  is  required  for  the  method.  Hot  pressing  results  in  good 
microstructure,  high  density,  moderate  grain  growth,  but  again  is  a  time- 
consuming  process. 

Similar  to  hot  pressing,  sinter  forging  samples  are  not  confined  laterally. 
High  strains  effectively  push  out  voids.  Although  a  hot  press  is  used,  in 
contrast  to  hot  pressing,  the  sample  in  sinter  forging  must  be  smaller  than 
the  die.  The  result  is  higher  density  and  minimal  grain  growth,  although 
the  process  usually  results  in  an  aligned  grain  structure  because  of  high 
strains. 

In  hot  isostatic  pressing,  inert  gas  pressures  up  to  30,000  psi  and 
temperatures  up  to  2000  °C  surround  the  SiC  compact.  A  hot  isostatic 
Press  (HIP)  unit  is  required  for  isostatic  pressing.  The  benefit  to  this  time- 
consuming  and  complicated  process  is  uniform  pressure  leading  to 
uniform  microstructure. 

Microwave  sintering  is  achieved  through  interaction  between  microwaves 
and  atoms  forming  the  material,  causing  rapid  heating  (reported  up  to 
1000  °C/min).  Through  a  combination  of  frequencies  from  2.5  to  85  GHz, 
densification  is  enhanced,  although  the  mechanisms  involved  are  unclear. 
Microwave  sintering  can  be  used  in  conjunction  with  pressure.  A 
microwave  furnace  is  needed.  The  process  produces  a  higher  density  than 
traditional  methods.  High  heating  rates  mean  shorter  sintering  times. 

Plasma  Assisted  Sintering  is  a  variant  of  microwave  sintering  in  which  the 
sample  is  surrounded  by  plasma  that  is  ignited  by  the  microwaves.  The 
plasma  promotes  enhanced  densification,  although  it  is  unknown  whether 
the  cause  is  the  plasma  or  the  high  temperature  inside  the  plasma. 

Spark  plasma  sintering  (SPS)  involves  passing  an  electrical  current  pulse 
though  the  sample  while  under  pressure  (30-50  MPa).  Spark  discharges  in 
the  voids  between  the  particles  are  assumed  to  cause  high  heating  rates 
(reported  up  to  600  °C/min)  leading  to  high  densification  rates.  This 
method  is  useful  for  nano-sized  particles.  This  method  is  used  because  of 
its  fast  sintering  times,  minimal  grain  growth,  and  high  density.  The  SPS 
sintering  method  using  an  arc  discharge  between  adjacent  particles 
requires  very  special  equipment.  Several  companies  make  these  systems 
including  FCT  in  Germany  (http://www.fct-svsteme.de/en/start/Home.html  and  SPS 
Syntax  in  Japan  (http://www.scm-sps.com/). 
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4.3  SiC  synthesis  experiments 

SiC  polycrystalline  samples  have  been  made  using  commercially  bought 
SiC  powder  and  sintering  additives.  The  powders  were  milled  to  fine  sizes, 
mixed  with  additives,  and  heated  to  high  temperatures.  By  the  end  of 
2010,  these  steps  were  being  optimized  to  form  a  dense  polycrystalline 
ceramic. 

The  as-received  powders  were  characterized  for  particle  size  and  density. 
After  planetary  milling  for  various  durations,  the  particle  size  was  again 
determined.  Representative  data  are  shown  in  Figure  20. 
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Figure  20.  Particle  size  for  various  milling  times. 

The  powders  were  then  placed  into  a  forming  mold  and  pressed  into  small 
disks,  which  were  subsequently  isostatically  cold-pressed  to  approximately 
45,000  psi.  The  samples  were  then  placed  in  graphite  crucibles  in  a 
vacuum  furnace  at  2000  °C  for  1  hour.  These  steps  are  shown  in  Figure  21. 

4.4  Product  characterization 

The  surface  and  interior  surface  were  characterized  by  SEM  images.  The 
samples  shown  in  Figure  22  are  not  fully  sintered.  They  do  not  have 
significant  density,  but  do  show  the  start  of  grain  diffusion.  Mechanical 
testing  has  yet  to  be  completed  for  any  samples. 
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a.  Composite  samples  after  dye  pressing,  before  CIP  (left)  and  within  a  protective  nitrile  glove 

ready  for  CIP  (right). 


b.  The  die  press  used  to  form  the 
discs  seen  in  figure  a. 


c.  The  CIP  unit  into  which  the 
discs  would  be  placed  after  being 
inserted  into  a  nitrile  glove  to 
prevent  liquid  infiltration. 


Figure  21.  a-b)  SiC  disks  before  CIP.  c)  Press  used  to  form  disks,  d)  CIP  unit. 


4.5  FY11  tasks 


The  tasks  of  acquiring  use  of  ceramic  synthesis  equipment,  ceramic 
synthesis  experiments,  and  making  refinements  to  the  simulation- 
synthesis  experimental  cycle  will  continue.  The  authors  hope  to  have 
access  to  an  SPS  unit  for  high-rate  sintering. 


Synthesis  methods  will  continue  to  be  explored,  particularly  attempting  to 
understand  what  occurs  at  each  of  the  many  processing  steps.  Powder 
synthesis  methods  including  ball  milling  and  mechanical-chemical  grind¬ 
ing  will  be  studied.  Powder  characterization  including  particle  size,  ag¬ 
glomeration,  density,  crystal  structure,  and  surface  chemistry  will  be 
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Figure  22.  Samples  after  sintering  step. 


investigated.  Experiments  for  mixing,  packing,  and  forming  of  powders 
will  be  performed.  The  sintering  techniques  chosen  will  be  employed;  and 
the  influences  of  particle  size,  particle  geometry,  grain  growth,  coarsening, 
pore  mobility,  green  body  density,  and  heating  schedules  on  the  resulting 
microstructure  will  be  documented. 

Using  the  knowledge  gained  on  the  methods  of  polycrystalline  ceramic 
synthesis,  SiC  samples  will  continue  be  made  and  refined.  The  objective  is 
to  make  a  polycrystalline  material  that  has  as  close  to  too  percent  density 
as  possible  and  contains  the  smallest  grains  possible.  This  work 
necessarily  includes  two  components:  synthesis  and  characterization.  Both 
components  are  straightforward,  but  take  considerable  time.  Processing, 
forming,  and  sintering  (the  three  steps  in  the  synthesis  part)  are  time 
consuming  and  require  attention  to  detail.  The  characterization 
component  includes  activities  such  as  scanning  electron  and  transmission 
electron  microscopy,  atomic  force  microscopy,  x-ray  diffraction,  density 
measurements,  and  mechanical  testing.  Each  of  these  activities  is  also  time 
consuming  and  will  need  to  be  repeated  for  each  sample  made. 

The  simulation-synthesis  experiment  cycle  will  begin  for  each  sample  by 
determining  the  microstructure  and  mechanical  properties,  which  will  be 
sent  to  the  modeling  and  simulation  teams,  who  will  model  the 
microstructure  for  simulation  of  mechanical  tests.  This  procedure  will 
allow  verification  that  the  algorithms  used  are  correct.  Once  verified,  the 
modeling  and  synthesis  teams  can  begin  to  work  together  on  the 
improvement  of  specific  properties  based  on  changes  in  modeled 
microstructure. 
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Once  an  appropriately  dense  SiC  sample  is  made,  silica-functionalized 
carbon  nanotubes  will  be  added  to  the  initial  powders.  A  description  of 
work  being  done  to  functionalize  CNTs  with  silica  is  given  later  in  this 
technical  report.  The  method  for  synthesizing  polycrystalline  ceramics  has 
been  chosen  to  be  powder  sintering.  The  sintering  methods  chosen  are 
liquid  phase  sintering,  hot  pressing,  and  SPS. 

4.6  Technical  gaps 

Fundamental  research  needs  to  be  conducted  to  determine  why  SPS 
works.  So  far  all  that  is  known  is  that  it  works  very  well  (i.e.,  produces 
dense  ceramics),  but  no  one  has  a  definitive  answer  as  to  what  occurs 
during  the  sintering  process.  The  authors  would  be  able  to  guide  their 
multiscale  modeling  more  accurately  if  they  knew  the  answer  to  this 
question,  particularly  the  physical  synthesis.  In  addition,  knowing  what 
happens  during  the  SPS  process  will  allow  better  prediction  of  what  might 
happen  to  CNTs  dispersed  in  the  matrix,  how  to  protect  them  from 
degradation,  how  to  allow  for  better  bonding,  and  how  to  make  full  use  of 
their  strength.  Ceramic  synthesis  and  sintering  have  been  around  for  quite 
some  time,  but  the  new  method  of  SPS  has  brought  about  the  capability  to 
make  dense  ceramics  with  small  grain  size.  How  these  small  grain  sizes 
interact  with  carbon  allotropes  is  one  question  that  needs  to  be 
investigated  during  this  project. 

4.7  Carbon  nanotube  functionalization 

A  review  of  literature  covering  the  functionalization  of  CNT  with  silica 
affords  several  different  methods,  each  with  their  own  drawbacks.  Li  et  al. 
(2009)  claim  that  their  process  does  not  damage  CNT  and  is  economical 
and  environmentally  friendly.  Their  procedure  relies  on  the  Stober  method 
for  producing  nanoscale  silica,  wherein  an  ammonia-activated 
cocondensation  reaction  occurs  between  ethanol,  hydroxyl  functional 
groups  on  CNTs,  and  tetraethyl  orthosilicate.  When  complete,  the  CNTs 
are  functionalized  with  a  silica  layer  2  to  10  nm  in  thickness  with  varying 
uniformity  (Figure  23).  The  thickness  and  uniformity  may  be  controlled  by 
various  reaction  parameters. 

First,  multi -walled  CNTs  or  single-walled  CNTs  are  functionalized  with 
hydroxyl  groups  50  -  250  mg  at  a  time  using  a  Harrick  Plasma  Cleaner 
(Figure  24).  The  instrument  provides  an  inductively  coupled  air  plasma  at 
about  150  mtorr  and  is  cycled  three  times  for  20  seconds  each  cycle.  These 
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CNTs  are  then  ultrasonically  dispersed  in  100  -  250  mL  of  ethanol  using  a 
Cole-Parmer  CPX130PB.  In  this  research,  sonication  imparts  a  total  of 
1.5  kJ  of  ultrasonic  energy  to  the  suspension  by  applying  30  percent  power 
for  5  min  using  a  V4  inch  probe.  Following  the  sonication  process,  4  - 
20  mL  of  28  percent  NH3  and  1-4  mL  of  tetraethyl  orthosilicate  are  add¬ 
ed  to  the  mixture,  which  is  then  covered,  stirred,  and  heated  at  65  °C  for 
0.5-4  hours.  The  mixture  is  then  immediately  filtered  over  a  membrane. 
To  obtain  a  suspension  of  the  silica-functionalized  CNT  in  water,  the  func¬ 
tionalized  product  may  be  added  to  water  and  sonicated. 


Figure  23.  Silica  functionalized  CNTs  made  at  ERDC  -CERL. 


Figure  24.  Harrick  plasma  cleaner  used  to 
functionalize  CNTs. 
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5  Modeling  the  Sintering  Process 

Sintering  is  a  process  by  which  powdered  SiC  is  coalesced  into  a  solid 
polycrystalline  mass  through  heating  and  sometimes  using  pressure.  The 
grain  size  distribution  of  the  polycrystalline  material  differs  from  that  of 
the  original  powder  depending  on  the  particular  sintering  process  and 
values  of  the  control  parameters  (see  Chapter  4).  In  modeling  the  sintering 
process,  two  questions  arise:  What  are  the  physical  conditions  required  to 
achieve  the  specified  material  properties,  and  can  those  conditions  be 
achieved  uniformly  throughout  the  body?  The  first  question  requires 
consideration  of  atomic-scale  and  mesoscale  processes.  The  second 
involves  modeling  the  energy  transfer  within  the  mass  using  continuum 
models. 

5.1  Particle-scale  view  of  sintering 

SiC  structural  materials  are  formed  by  either  solid-  or  liquid-state 
sintering  (Niu  et  al.  2011  and  Biswas  2009).  Solid-phase  sintering 
includes  two  stages:  necks  form  at  the  contact  points  between  particles, 
then  the  grain  boundary  diffusion  causes  shrinkage,  and  the  material  is 
densified.  Liquid-phase  sintering  is  the  process  of  including  an  additive  in 
the  powder  that  will  melt  before  the  matrix  phase.  The  process  of  liquid 
phase  sintering  has  three  stages:  rearrangement,  solution-precipitation, 
and  final  densifi cation.  Mesoscale  models  have  been  developed  to  study 
the  evolution  of  microstructure  during  both  solid  sintering  (Braginsky  et 
al.  2005)  and  liquid  sintering  (Lee  et  al.  2007). 

5.1.1  Mesocale  modeling  of  solid-state  sintering 

Braginsky  et  al.  (2005)  used  a  kinetic  Monte  Carlo  (KMC)  model  to 
simulate  two-dimensional  microstructure  evolution  during  solid-state 
sintering.  The  model  can  capture  grain  growth  by  short-range  diffusion 
from  one  side  of  a  grain  boundary  to  the  other,  and  it  can  also  model  long- 
range  diffusion  of  pores  by  surface  diffusion  and  vacancy  annihilation  at 
grain  boundaries. 

In  the  model  by  Braginsky  et  al.  (2005),  grains,  pores  and  vacancies  are 
represented  by  an  assemblage  of  points.  Each  point  is  assigned  a  state 
corresponding  to  a  unique  integer  value,  qgram  =  [1,2, ...Q] ,  that 
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corresponds  to  Q  different  crystallographic  orientations.  Each  pore  site  is 
assigned  only  one  state  qp0re  =-i.  Each  grain  consists  of  contiguous  sites 
with  the  same  orientation  (i.e.,  same  integer  value).  Grain  boundaries  are 
therefore  represented  as  the  surfaces  between  sites  of  different  grains,  and 
grain  surfaces  are  the  boundary  between  grain  sites  and  pore  sites.  A 
vacancy  is  defined  in  this  model  as  a  single  pore  site  that  is  entirely 
surrounded  by  grains. 

The  primary  driving  force  for  this  model  may  be  expressed  as  the 
reduction  in  interfacial  free  energy.  The  equation  of  state  is  thus  the 
summation  of  interaction  energy  of  the  system,  defined  as 

E  =  YiU^UhJ[^-^qi,qj-)\  (8) 


where 


N 

n 

Jij 

<P 

Qj 

Siquqj) 


total  number  of  points 

number  of  adjacent  neighbor  points 

neighbor  interaction  energy  between  sites  i  and  j 

state  of  the  grain  or  pore  at  site  i 

state  of  the  nearest  neighbor  at  site/ 

Kronecker  delta  with  8 (gz  =  qj)  =  l  and  8 (q,  *  qj)  =  o 


The  standard  metropolis  algorithm  is  then  used  to  perform  the  grain 
growth  step  with  probability  P  given  by 


exp 


AE 


kBT 


for  AE  >  0 


1  for  AE  <  0 


(9) 


where 

E  =  total  energy  using  Equations  8  and  g 
Kb  =  Boltzmann  constant 
T  =  temperature 

Figure  25  provides  an  example  utilizing  an  in-house-developed  2D  KMC 
model  involving  the  microstructural  evolution  of  a  sintered  powder 
compact  over  a  total  duration  of  1  million  Monte  Carlo  steps  (MCS).  The 
time  in  a  KMC  model  is  measured  in  MCS,  such  that  1  MCS  corresponds  to 
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N  attempted  grain  growth  steps.  As  indicated,  the  initial  microstructure 
shows  a  large  numbers  of  fine  grains  and  pores.  As  the  simulation 
continues,  the  number  of  pores  decreases  and  the  grain  sizes  are  several 
times  larger  than  those  at  the  initial  stages  of  sintering. 


a.  t=0  MCS  b.  t=1000  M  c.  t=l  E6  MCS 

Figure  25.  Microstructure  evolution  of  a  random  packing  nonuniform  powder. 

5.1.2  Mesoscale  modeling  of  liquid  phase  sintering 

Lee  et  al.  (2007)  performed  three-dimensional  simulations  of  isotropic 
coarsening  in  liquid  phase  sintering.  In  liquid  phase  sintering,  particle 
coarsening  occurs  by  mass  transport  through  the  liquid  surrounding  solid 
particle.  The  driving  force  is  determined  by  the  solubility  of  particles  (Lee 
et  al.  2007): 


In 


2  My 
RgTpr 


(10) 


where 

Sco  =  solubility  of  a  solid  particle  with  an  infinite  radius 
M  =  molecular  weight  of  the  solid  with  density  p 
y  =  interfacial  free  energy  between  the  solid  and  liquid  phases 
(assumed  to  be  isotropic) 

Rg  =  gas  constant 
p  =  density  of  solid 
r  =  radius  of  solid  particle 
T  =  absolute  temperature 

Monte  Carlo  simulations  also  are  used  to  model  the  grain  growth  in  the 
liquid-phase  sintering.  Similar  to  the  model  used  for  solid-phase  sintering, 
an  assembly  of  voxels  populates  the  solid  phase  and  liquid  phase.  Each 
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solid  voxel  is  assigned  a  spin  number  based  on  the  local  particle 
orientation,  and  each  liquid  voxel  is  assigned  a  spin  number  of  -l.  The 
kinetic  Monte  Carlo  simulations  are  performed  to  minimize  the  interface 
energy. 

5.1.3  Technical  gaps 

Both  solid-phase  sintering  and  liquid-phase  sintering  are  complex 
processes.  Many  factors  such  as  temperature,  pressure,  and  sintering 
additives  affect  the  development  of  microstructure  (Maizza  et  al.  2007; 
Dillon  and  Harmer  2008;  Pan  and  Gong  2010).  Current  mesoscale 
sintering  models  typically  address  only  grain  growth,  pore  migration,  and 
vacancy  annihilation  at  grain  boundaries  by  minimizing  the  total  surface 
energy.  This  will  often  include  grain  boundary  mobility  effects,  which  take 
into  account  the  curvature-driven  grain  growth  of  boundaries  toward  their 
centers  of  curvature.  The  effects  of  temperature,  pressure,  and  additives 
are  not  included  in  the  current  models.  The  size,  shape,  and  distribution  of 
initial  powders  also  affect  the  microstructure  development.  The  powders 
are  usually  compacted  into  green  bodies  before  the  high  temperature  is 
applied.  The  current  models  usually  use  spherical  particles  with  random 
radius  and  packing  to  represent  powders,  and  do  not  include  the 
compaction  process.  Additives  can  affect  grain  growth  and  phase 
transformation  (Pan  and  Gong  2010).  The  effects  of  additives  are  not 
included  in  the  grain  growth  models. 

5.2  Modeling  the  sintering  environment 

The  continuum  scale  modeling  of  SiC  synthesis  due  to  heat  transfer  is 
complex.  This  complexity  is  due  largely  to  the  varied  and  closely  coupled 
heat  transfer  mechanisms  that  occur  during  the  heating  process.  These 
mechanisms  include  conduction,  convection,  and  radiation.  Additionally, 
the  primary  heat  sources  are  electromagnetic  in  nature,  and  as  such, 
require  the  numerical  solution  of  Maxwell’s  equations. 

Because  of  this  modeling  complexity,  a  “bottom-up”  approach  was  selected 
for  this  work.  The  approach,  as  its  name  suggests,  begins  with  the 
successful  development,  execution,  and  validation  of  the  simplest  models 
and  continues,  iteratively,  toward  models  of  higher  complexity.  The  first  of 
these  models  consists  of  an  infinite  rod  of  constant  radius  and  current 
density,  for  which  an  analytical  solution  exists.  This  model  incorporates  a 
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simplified  version  of  the  electromagnetic  equations  and  is  based  on  a 
solution  to  the  electric  scalar  and  vector  potentials. 

Following  the  infinite  rod  problem,  the  modified  Lely  method  was  selected 
as  the  next  phase  in  the  development  process.  The  method  involves 
sublimation  of  SiC  from  a  hot  powder  source,  transport  of  vapor  through 
an  inert  gas  environment,  and  condensation  on  a  seed  that  is  colder  than 
the  source.  The  deposition  continues  until  a  bulk  crystal  of  reasonable  size 
is  grown.  The  method,  in  addition  to  low-frequency  induction  heating, 
requires  the  numerical  modeling  of  several  heat  transfer  mechanisms, 
including  conduction,  radiation,  and  potentially  convection.  The 
numerical  modeling  effort  is  further  complicated  by  the  fact  that  the 
growth  apparatus  is  nonhomogenous,  but  is  composed  of  several  different 
materials  types,  with  temperature-dependent  properties. 

The  final  phase  of  the  bottom-up  approach  is  still  under  initial 
development.  This  phase  attempts  to  predict  the  temperature  and  stress 
field  evolution  during  the  SPS  process.  SPS  is  a  powder  consolidation 
technique  that  employs  pulsed  electric  current  and  mechanical  pressure  to 
assist  in  the  sintering  process  (Munoz  2010).  This  coupled  electrical- 
thermal-mechanical  analysis  will  not  only  incorporate  all  of  the 
aforementioned  modeling  capabilities  encountered  in  both  the  infinite  rod 
problem  and  modified  Lely  method,  but  will  also  require  modeling  the 
stress  distribution  due  to  mechanical  inputs. 

The  primary  objective  of  this  research  is  thus  to  provide  a  status  of 
capabilities  toward  the  ultimate  goal  of  modeling  the  continuum  SPS 
synthesis  process.  From  an  overall  context,  this  includes  providing 
sufficiently  accurate  thermal  and  stress  field  calculations  to  be  used  by 
other  members  of  the  Advanced  Material  Initiative  (AMI)  team. 

5.2.1  Infinite  rod 

Process.  As  a  first  approximation  to  the  solution  of  the  set  of  electromag¬ 
netic  equations,  a  static  electromagnetic  problem,  with  a  known  analytical 
solution  was  selected  (Sass-Tisovskaya,  M.  2009).  Taking  into  account  the 
assumptions  of  electroneutrality,  quasi-steady  state,  and  neglecting  the 
Hall  and  induction  currents,  Maxwell’s  electromagnetic  equations  take  the 
following  form: 


V-j  =  0 


(11) 
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VB  =  0 

(12) 

V-E  =  0 

(13) 

V  x  B  =  fuj 

(14) 

V  xE  =  0 

(15) 

j  =  ^cE  , 

(16) 

where 

j  =  current  density 
B  =  magnetic  flux  density 
E  =  electric  field 

|i0  =  magnetic  permeability  in  a  vacuum 
ac  =  electric  conductivity. 

Combining  the  set  of  Equations  12-17,  the  potential  form  of  the 
electromagnetic  field  may  be  written  as: 


V-(crV^)  =  0 

V2A  =  acEoV</>  ’  (17) 

where  cp  and  A  are  the  electric  scalar  potential  and  magnetic  vector 
potential,  respectively. 

Consider  an  infinite  rod  of  constant  radius  r0,  and  a  constant  current 
density  running  parallel  with  the  rod  axis.  The  electrical  conductivity  is 
constant  inside  the  rod,  representing  a  hot  argon  gas  column  at  a  fixed 
temperature  of  T  =  10,600  K. 

Inside  the  cylinder,  the  analytic  solution  of  the  magnetic  vector  potential 
(Ax  in)  takes  the  form 


Ay,  in  A) 


EoJA 


■  if  r<r0 


4 


(18) 
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where  A0  is  the  magnetic  vector  potential  at  the  fixed  temperature. 

Outside  the  cylinder,  the  magnetic  vector  potential  (Ax  out)  takes  the  form 

A,ou,  =Ao~  [0-5  +  In (r  /  r0 )] ,  if  r  >  r0,  (19) 

where  jx  can  be  found  via  the  current  intensity  I  from: 

j  =  -b  (20) 

7Tr0 

Because  of  the  symmetry  of  the  problem,  the  geometry  was  modeled  as  a 
two-dimensional,  axisymmetric  wedge,  and  a  uniform  hexagonal  mesh 
was  used.  The  boundary  conditions  for  the  magnetic  scalar  potential 
included  zero  gradients  normal  to  all  domain  surfaces  except  the  outer 
surface  of  the  wedge.  The  electric  scalar  potential  cp  was  prescribed 
similarly  except  for  the  two  ends  of  the  rod,  which  were  fixed  at  cp  =  oV 
and  cp  =  707V. 

Results.  The  numerical  results  are  presented  in  Figures  26  and  27.  The 
isolines  corresponding  to  the  electric  potential  (see  Figure  26)  run  parallel 
to  each  other  inside  the  rod  and  curve  symmetrically  outward  from  the 
rod.  This  pattern  is  due  primarily  to  the  prescribed  boundary  conditions 
for  the  electric  potential.  Figure  27  compares  the  numerical  and  the 
analytic  solutions  for  the  x-component  of  the  magnetic  vector  potential 
along  the  radius  of  the  domain  (x  =  0.5E-2  m)  and  shows  good  agreement 
with  the  analytical  results  computed  from  Equations  18  and  19. 

5.2.2  Modified  Lely  method 

Sublimation  growth  via  the  modified  Lely  method  has  been  one  of  the 
most  successful  and  most  widely  used  SiC  single  crystal  growth  techniques 
of  recent  years  (Klein  et  al.  2000).  The  process  involves  sublimation  of  SiC 
from  a  hot  powder  source,  transport  of  vapor  through  an  inert  gas 
environment,  and  condensation  on  a  seed  that  is  colder  than  the  source. 
The  deposition  continues  until  a  bulk  crystal  of  reasonable  size  is  grown. 

As  shown  in  Figure  28,  an  SiC  growth  apparatus  typically  consists  of  a 
Radio  Frequency  (RF)  copper  coil,  a  quartz  tube,  a  graphite  susceptor, 
insulation,  and  an  RF  induction-heated  graphite  crucible  containing 
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Figure  26.  Electric  potential  isolines  together  with  current 
density  vectors. 


Figure  27.  Numerical  and  analytical  results. 
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Figure  28.  Schematic  of  the  RF  induction  heating  system  for  growth  of  SiC  crystals. 


polycrystalline  SiC  source  powder  and  a  single  crystalline  SiC  seed.  The 
SiC  source  powder  is  heated  using  RF  induction  heating  via  copper  coils. 
The  heating  is  generated  from  a  current  of  approximately  10  kHz  and  1200 
amps.  Due  to  the  high  electrical  conductivity  of  the  graphite  susceptor,  the 
time-harmonic  electromagnetic  field  induces  eddy  currents  therein,  and 
heat  is  generated  via  the  Joule  effect.  In  general,  the  electromagnetic  wave 
penetrates  the  susceptor  within  a  skin  depth  of  only  a  few  millimeters,  but 
is  generally  adequate  to  allow  for  sufficient  heat  generation  within  the 
inner  growth  region  via  conduction  and  radiation.  Convection  may  also 
contribute  in  larger  growth  systems  (Chen  et  al.  2001).  The  SiC  charge 
powder  is  placed  within  the  lower,  heated  portion  of  the  crucible,  and  the 
seed  crystal  is  placed  above  and  is  cooled  with  a  blind  hole.  This 
temperature  difference  allows  sublimation  to  occur  at  the  source  and 
crystallization  to  form  at  the  seed. 

The  growth  rate  of  the  crystal  has  been  shown  in  numerous  studies  to  be 
strongly  linked  to  the  temperature  distribution,  particularly  with  respect  to 
the  temperature  difference  between  the  source  and  seed  (Barrett  et  al. 
1993;  Kaneko  1993;  Hofman  et  al.  1995;  Buhner  et  al.  1999;  and  Ma  et  al. 
2000;).  Since  SiC  dissociates  completely  into  liquid  silicon  and  solid 
carbon  at  3150  K,  the  growth  temperature  must  remain  below  this  upper 
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limit  (Lilov  1993).  The  temperature  distribution  inside  the  growth  system 
is  determined  via  complex  interactions  among  the  various  heat  transfer 
mechanisms,  and  occurs  within  materials  of  different  phases,  including 
gases,  solids,  and  porous  media.  The  temperature  further  depends  on 
several  configuration  control  parameters,  including  coil  position,  heating 
voltage,  and  pressure. 

Since  favorable  growth  conditions  are  critically  linked  to  the  temperature 
distribution,  its  optimization  is  an  important  goal  of  apparatus  design  and 
control.  However,  because  of  the  high  temperatures  (approximately 
3000  K),  experimental  verification  of  the  correlation  between  the 
properties  of  the  crystal,  the  temperature  distribution,  and  external 
control  can  be  difficult.  As  a  result,  accurate  numerical  simulations  have 
become  highly  desirable. 


The  primary  objective  of  this  phase  of  the  work  is  to  develop  a  coupled, 
high-resolution  heat  transport  model  to  predict  the  electromagnetic  field 
and  temperature  variation  due  to  RF  heating  during  SiC  crystal  synthesis 
using  the  modified  Lely  method.  The  model  proposed  here  can  be  used  to 
investigate  the  effects  of  variable  current  intensity  and  coil  position  on  the 
temperature  distribution  within  the  growth  chamber. 


5.2.3  Electromagnetic  field 

The  electromagnetic  field  produced  by  RF  induction  heating  can  be 
calculated  using  Maxwell’s  equations.  For  low-frequency  if  <  1  MHz) 
induction  heating,  Maxwell’s  equations  can  be  simplified  using  the  quasi- 
steady-state  approximations.  Assuming  the  current  in  the  coil  is  time- 
harmonic,  the  magnetic  flux  density  can  be  expressed  as  the  curl  of  a 
magnetic  vector  potential,  B  =  V  x  A.  The  Maxwell  equations  can  then  be 
written  in  terms  of  the  vector  potential,  A  (Jackson  1998): 


Vx 


— VxA 

vA Lm 


d2A  <9A  . 

+e”'aF+^m-1 


coil 


(21) 


where 

/Lim  =  magnetic  permeability 
A  =  magnetic  potential  vector 
e  m  =  permittivity 
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oc  =  electrical  conductivity 
J coil  =  current  density  in  the  coil. 

If  the  coil  and  the  electromagnetic  field  are  assumed  axisymmetric,  both 
the  magnetic  potential  vector  A  and  the  current  density  Jcoiz  will  have  only 
one  angular  component  with  an  exponential  form,  such  that 

A=U,e“+°^'“l;  -u={vi“+°v-‘“}  (22) 

where  i  is  the  complex  unit,  and  co  is  the  angular  frequency.  The  final 
equation  for  vector  potential,  A,  is  obtained  by  substituting  Equation  23 
into  22: 
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(23) 


The  second  term  in  the  left  side  of  Equation  24  is  negligible.  The  following 
boundary  conditions  can  be  used  to  solve  this  equation, 

A0  =0atr  =  0and(r2  +z2)-^oo  (24) 

which  is  based  on  the  conditions  that  magnetic  flux  density  B  is  axisym¬ 
metric  and  the  computational  domain  is  sufficiently  large  (Wu  et  al.  2005). 
After  Equation  24  is  solved  for  the  vector  potential  A,  the  generated  heat 
power  in  the  graphite  susceptor  can  be  obtained  using  the  principles  of 
eddy  current: 


Qeddy  =|^c«2Ur+A2) 


(25) 


5.2.4  Radiation  heat  transfer 

Diffuse  gray  radiation  is  modeled  using  heat  flux  boundary  conditions  at 
the  relevant  heat  exchange  surfaces.  For  this  work,  these  include  the 
surfaces  pertaining  to  the  SiC  source  powder  and  the  growth  chamber 
filled  with  inert  gas.  Additionally,  radiation  heat  transfer  modeling  may 
also  be  applied  to  the  top  and  bottom  cooling  holes. 
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In  general,  heat  flux  can  consist  of  idealized  radiation: 

=“(*- -C) 


(26) 


where 

o  =  Stefan-Boltzmann  constant 
n  =  surface  normal 

8  =  surface  emissivity  (considered  herein  constant  for  each 
surface) 

Text  =  external  temperature. 

If  the  surface  k  is  receiving  radiation  from  another  surface  i  in  the  system, 
then  the  heat  flux  may  be  represented  using  Elmer  software 
(http://www.csc.fi/english/pages/elmer)  as 


OE,. 
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\£k  j=l 


EG*£iTiAi 


(27) 


where  the  parameters  A,  andA/t  refer  to  the./  and  k  specific  surface  areas, 
respectively,  and  the  factors  Gik  represent  the  Gebhardt  factor: 


Q* 


jk 


Gjl<  £  :A  ;<jT? 


(28) 


Here  Qjk  is  the  heat  transferred  from  surface  j  to  k,  and  is  computed  from 
the  product  of  the  view  factors  Fjk  and  the  radiosity  J,.- 


p  -  9a 

jk  A.J  . 


(29) 


5.2.5  Material  properties 

The  temperature-dependent  material  properties  were  taken  from  the  liter¬ 
ature  (Chen  et  al.  2001).  The  SiC  source  material  was  modeled  with  an  av¬ 
erage  grain  size  dp  of  125  pm  (Sasaki  et  al.  1998)  and  was  treated  as  a  po¬ 
rous  medium  with  porosity  &p  and  effective  specific  heat  (Cp) 
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(pcp)eff=(l-sp)(pcp)sic+£p(pcp)gas  (30) 

where  p  is  the  density,  and  a  porosity  of  0.4  was  utilized  (Chen  et  al. 

2001). 

The  heat  transfer  inside  the  charge  consists  of  both  radiative  and 
conductive  components;  as  such,  an  expression  for  the  effective  thermal 
conductivity  must  account  for  both  of  these  components.  From  Kansa  et 
al.  (1977), 


keff=^-SP^kSiC 


-a4  oT3dD 
3  p 


(31) 


5.2.6  Convective  heat  transfer 

The  relative  importance  of  convective  heat  transfer  (particularly  natural 
convection  due  to  buoyancy  effects),  within  the  context  of  the  aforemen¬ 
tioned  crucible  dimensions,  as  well  as  the  processing  temperature  and 
pressure  conditions,  may  be  approximated  by  the  dimensionless  Rayleigh 
number: 

Ra  =  Gr  *  Pr  =  —  AT!D3  (32) 

va 

where 

Gr  =  Grashoff  number 
Pr  =  Prandtle  number 
P  =  thermal  expansion  coefficient 
v  =  kinematic  viscosity 
a  =  thermal  diffusivity 
g  =  acceleration  due  to  gravity 

AT  =  temperature  difference  between  the  crucible  wall  and  the 
quiescent  gas 

Dc  =  diameter  of  the  crucible. 

For  a  typical  growth  system  composed  of  argon  gas,  at  a  pressure  of  3.0E4 
Pa,  temperature  of  2900  K,  and  temperature  difference  of  30  K  (Chen  et 
al.  2001),  the  Rayleigh  number  may  be  anywhere  between  1.98  and  15.84 
for  representative  (industry  standard)  crucible  diameters  between  50  and 
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100  mm,  respectively.  Since  this  is  significantly  less  than  the  critical 
Raleigh  number  Rac,  it  is  evident  that  the  heat  transfer  by  natural 
convection  can  be  neglected  in  these  studies  (Reid  and  Harris  1958).  If 
however,  the  crucible  diameter  is  increased  (i.e.,  by  a  factor  >  4),  the 
relative  importance  of  the  buoyancy  effects  requires  the  inclusion  of 
convective  transport  modeling. 

5.2.7  Numerical  method 

In  concert  with  Equations  26,  32,  and  33,  the  equation  governing  the 
evolution  of  temperature  for  the  growth  system  can  be  stated  as 

=  +  (33) 

Additionally,  the  boundary  conditions  corresponding  to  surface  heat  flux 
(Equation  28)  and  magnetic  vector  potential  (Equation  25)  are  applicable. 

The  variational  form  of  Equation  35,  after  integration  by  parts  of  the 
conduction  term  and  the  application  of  the  divergence  theorem  yields 

f  — f  keff  vr-Vu.-dfl 
Jn  dt  {  (CpPlff 


=  f  -^-Vidn+  f 


kpPlff 


keffVT)rw,dT 


(34) 


where 

V;  =  basis  function 

H,T  =  elemental  volume  and  its  enclosing  surface,  respectively 
n  =  unit  normal  vector  to  the  surface. 


The  linearized  partial  differential  equation  for  the  temperature  (T)  in  the 
variational  formulation  is  stated  as: 


dT. 

M, — -  +  A.T.  =  F- 


dt 


V  J 


(35) 


where 
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Mij  =  mass  matrix 
Aij  =  stiffness  matrix 
Fi  =  force  vector. 

Comparing  Equation  34  and  35  leads  to  the  following: 


Mii  =  f  V,Vidn 


a 


W 


(cDp)e 


-Vu.  •  Vu.dff 


f  Qeddy  Vv  dn  +  r  Qn_  dr 
j  (CpPlff  r  CpP 


(36) 


where  qn  =  keffVT  n  is  the  heat  flux  perpendicular  to  the  surface  normal. 


The  numerical  integration  of  Equation  35  over  each  quadrilateral  element 
is  approximated  by  the  sum  over  all  Gauss-point  contributions,  and  solved 
using  the  open-source,  multiphysics  simulation  software  Elmer 
f http://www.csc.fi/engiish/pages/eimer).  Time-stepping  was  performed  using  a 
second-order  Crank  Nicolson  scheme  and  was  solved  to  steady-state 
convergence  with  an  error  tolerance  of  approximately  1.0  E-5. 


5.2.8  Results 

A  computational  domain  of  0.25  by  0.4  m  was  selected  for  this  work,  and 
its  size  verified  to  allow  for  a  sufficiently  diminishing  magnetic  potential  at 
the  far-field  boundaries,  where,  in  accordance  with  Equation  25,  a  zero 
magnetic  potential  condition  was  imposed.  A  grid  composed  of  20,916 
unstructured  elements  was  used,  and  clustering  was  performed  around  the 
areas  of  high-temperature  gradients  (Figure  29).  Figure  30  shows  the  in- 
phase  ( Ar )  and  out-of-phase  (A,)  contours  of  magnetic  potential,  for  five 
turns  of  coil,  a  current  of  1200A,  and  a  frequency  of  10  kHz.  As  indicated, 
the  contours  are  concentrated  along  the  outer  portion  of  the  graphite 
susceptor  and  tend  to  bend  inwards  along  the  top  and  bottom  of  the 
susceptor.  The  in-phase  component  reaches  its  maximum  strength 
(-0.0018  Webers/m)  around  the  coil,  with  its  shape  being  highly  influ¬ 
enced  by  the  geometry  of  the  coil.  A,-,  in  contrast,  shows  a  maximum  mag¬ 
nitude  of  approximately  0.00017  Webers/m  and  is  located  both  near  the 
coil  and  within  the  graphite  susceptor.  Of  particular  interest  is  the  fact  that 
both  A;  andAr  potentials  show  comparable  strength  within  the  graphite 
susceptor  and  suggest  the  need  to  include  both  components. 
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Figure  29.  Computational  grid  for  finite  element  method. 
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Figure  30.  Contours  of  magnetic  potential,  out  of  phase  (Ai)  and  in-phase  (Ai), 
for  a  current  of  1200  Amps  and  a  frequency  of  10  kFIz. 
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Figure  31  compares  eddy  current  heat  generation  with  the  existing 
literature  (Chen  et  al.  2001).  Specifically,  the  heat  generation  within  the 
SiC  powder  charge,  the  crucible  (modeled  as  an  inert  gas),  and  the 
graphite  susceptor  are  shown  along  the  radial  direction  (r)  at  the  various 
heights  (z),  corresponding  to  z  =  0.0,  z  =  1.3 Rs,  and  z  =  2.6 Rs  (where  Rs  is 
the  outside  radius  of  the  susceptor,  Rs  =  0.07  m).  As  indicated,  the  heat 
power  is  generated  primarily  in  the  outer  portion  of  the  graphite  susceptor 
and  reaches  an  approximate  value  of  8.7  E6  Joules/m3.  The  discontinuity 
between  r/Rs  =  0.6  and  r/Rs  -0.8  occurs  as  a  result  of  the  intermediate 
inert  gas  material.  The  heat  generated  within  the  SiC  charge  at  z  =  i.3-Rs  is 
somewhat  negligible  compared  with  that  generated  in  the  graphite 
susceptor  at  z=o  and  2.6 Rs,  showing  an  approximate  difference  of  two 
orders  of  magnitude.  Overall,  the  numerical  results  of  this  study  compare 
well  with  the  results  of  Chen  et  al.  (2001). 


Figure  31.  Comparison  of  generated  heat  power  along  the  radial  direction  at  different  heights 

with  results  from  Chen  et  al.  (2000). 


Figure  32  shows  the  temperature  distribution  for  all  components  minus 
the  coil  apparatus.  As  indicated,  the  maximum  temperature  (-3075  K)  ex¬ 
ists  within  the  graphite  susceptor,  along  the  geometric  center  of  the  induc¬ 
tion  coil,  and  extends  radially  inward  towards  the  SiC  source  powder.  The 
coil  geometry  is  thus  responsible  for  the  location  of  the  highest  tempera- 
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ture  distributions  within  the  source  powder,  and  therefore  the  first  tem¬ 
perature  distribution  has  been  documented  in  several  studies  (e.g.,  Chen  et 
al.  2001),  which  show  these  attributes  to  be  highly  correlated. 


Figure  32.  Temperature  contours  for  a 
current  of  1200  A  and  frequency  of 
10  kHz. 


The  process  of  sublimation  continues  as  a  result  of  the  temperature 
difference  that  exists  between  the  source  powder  and  the  seed.  As  shown 
in  Figure  32,  the  temperature  at  the  top  of  the  growth  chamber  {Ttop)  is 
cooled  by  a  hole  exiting  to  the  outside  air.  This  hole  helps  maintain  the 
positive  temperature  difference  between  the  source  and  the  seed  and  thus 
facilitates  the  continued  sublimation  growth  at  the  seed.  This  temperature 
difference  is  approximately  80  K  and  allows  for  a  sufficiently  high  growth 
rate.  Further,  this  hole  helps  to  maintain  a  proper  radial  temperature 
gradient  at  the  seed  (i.e.,  thermal  variations  of  less  than  10K  at  the  seed 
surface)  that  ensures  low  thermal  stresses,  which  subsequently  assist  in 
minimizing  dislocations  and  micropipes  (Chen  et  al.  2001). 

The  effect  of  radiation  modeling  on  the  temperature  distribution  can  be 
seen  in  Figure  33.  The  temperature  difference  between  neighboring 
isolines  is  10K.  As  shown,  neglecting  radiation-modeling  results  in  large 
heat  gradients  inside  the  growth  chamber,  where  experiments  have  shown 
the  temperature  distribution  to  be  nearly  homogeneous.  As  stated 
previously,  these  erroneous  predictions  can  lead  to  false  conclusions  not 
only  with  respect  to  the  formation  of  dislocations  and  micropipes  resulting 
from  uncharacteristically  high  thermal  stresses  but  also  with  respect  to  SiC 
crystal  growth  rates  (Chen  et  al.  2001). 
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Figure  33.  Effect  of  applying  radiation  to  the  numerical  simulation. 


5.3  Conclusions 

A  comprehensive  process  model  for  SiC  growth  by  the  modified  Lely 
method  has  been  developed  that  accounts  for  axisymmetric  geometries, 
induction  heating,  radiation,  and  conductive  heat  transfer.  The  model  has 
been  successfully  applied  to  the  simulation  of  SiC  growth  and  compared 
with  existing  numerical  experiments  (Chen  et  al.  2001).  Computations  of 
the  magnetic  potential  showed  the  relative  importance  of  including  both 
the  in-phase  and  out-of-phase  components,  with  comparable  magnitudes 
existing  within  the  graphite  susceptor.  The  graphite  susceptor  with  its 
relatively  high  conductivity  was  shown  to  function  as  a  barrier  to  the 
electromagnetic  field,  allowing  only  a  negligible  amount  of  eddy-induced 
heating  within  the  inner  crucible  and  growth  region.  By  contrast,  a 
significant  amount  of  eddy  current  energy  was  created  in  the  outer  portion 
of  the  susceptor  within  a  relatively  small  skin  depth.  The  maximum 
temperature  distribution  was  found  to  exist  in  the  graphite  susceptor,  and 
at  the  level  of  the  geometric  center  of  the  induction  coils.  Within  the  SiC 
charge,  the  temperature  was  highest  in  the  central  region  near  the  crucible 
wall.  The  positive  temperature  difference  between  the  SiC  charge  and  the 
ambient  air-cooled  seed  allows  for  the  sublimation  in  the  charge  and 
deposition  on  the  seed.  The  effect  of  radiation  modeling  was  significant, 
particularly  given  the  high-temperature  ranges  appropriate  to  this  study, 
and  resulted  in  a  substantial  decrease  in  temperature  gradients  within  the 
growth  chamber. 
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Future  simulations  will  evaluate  the  effects  of  buoyancy  by  convection 
caused  by  increasing  the  growth  chamber  dimensions.  These  simulations 
will  necessarily  require  the  coupling  of  the  current  model  with  the 
compressible  Navier-Stokes  equations. 

5.3.1  Spark  plasma  sintering 

This  phase  of  the  work  attempts  to  predict  the  temperature  and  stress  field 
evolution  during  the  Spark  Plasma  Sintering  (SPS)  process.  SPS  is  a 
powder  consolidation  technique  that  employs  pulsed  electric  current  and 
mechanical  pressure  to  assist  in  the  sintering  process  (Munoz  2010). 

While  several  recent  studies  have  used  the  finite  element  method  to 
investigate  the  current  and  temperature  distributions  in  the  SPS  process 
(Anselmi-Tamburini  et  al.  2005;  Antou  et  al.  2009;  McWilliams  and 
Zavaliangos  2008;  Vanmeensel  et  al.  2005;  Wang  et  al.  2007;  and 
Zavaliangos  et  al.  2004),  the  mechanical  input  and  the  resulting  stress 
distributions  have  largely  been  neglected.  Experimental  studies  show  this 
to  be  a  significant  omission  because  the  processing  pressure  is  a  crucially 
important  densification  mechanism  (Anselmi-Tamburini  et  al.  2006); 
(Munir  et  al.  2006).  Only  a  few  exceptions  exist  that  include  coupled 
electrical-thermal-mechanical  analysis  (e.g.,  Munoz  et  al.  2010;  Wang  et 
al.  2007;  and  Antou  et  al.  2009). 

This  work  will  also  incorporate  a  coupled  electrical -thermal-mechanical 
analysis  and  will  use  one  or  more  of  the  references  cited  in  the  preceding 
paragraphs  for  validation.  Clearly  the  modeling  effort  will  not  only 
incorporate  the  capabilities  demonstrated  in  both  the  infinite  rod  problem 
and  modified  Lely  method,  but  will  also  require  the  need  for  modeling  the 
stress  distribution  due  to  a  mechanical  input  and  will  thus  require 
solutions  to  the  strain  displacement  relations. 

The  outputs  that  will  be  provided  to  other  AMI  researchers  from  the  con¬ 
tinuum  SPS  model  will  therefore  include  the  thermal  and  stress  fields  con¬ 
tained  within  the  growth  chamber. 


ERDC  TR-12-9 


68 


6  Material  Diagnostics 

A  large  amount  of  research  has  been  performed  in  the  scientific 
community  to  understand  the  mechanical  properties  of  various  materials, 
and  how  they  may  be  influenced  by  factors  such  as  stress  states,  length 
scales,  and  the  internal  structure  of  the  material.  As  discussed  in 
Horstemeyer  (2009)  on  multiscale  modeling,  Bazant  and  Chen  (1997) 
summarized  the  research  to  date  on  length  scale  effects  on  properties.  As 
early  as  the  1500’s  DaVinci  claimed  two  ropes  of  equal  thickness  and 
material  but  with  different  lengths  exhibited  different  tensile  strengths: 
the  longer  of  the  two  ropes  would  not  be  as  strong  as  the  shorter  rope.  As 
for  compressive  strength,  in  the  1700’s,  Euler  examined  column  length 
with  respect  to  buckling.  Research  performed  by  Cauchy  in  the  1800’s  and 
Bridgman  (1923)  in  the  1900’s  examined  stress  state  effects  on  materials. 
More  currently,  as  technology  progressed  and  diagnostic  capabilities 
improved,  Griffith  (1921)  pioneered  fracture  mechanics  developing 
Griffith’s  criterion,  which  is  most  applicable  to  brittle  materials,  while 
Irwin’s  (1957)  research  group  included  plasticity  to  better  predict  the 
fracture  behavior  of  ductile  materials.  Additionally,  Olson  (1998,  2000) 
identified  how  various  structures  such  as  grains,  defects,  and  inclusions  in 
a  certain  material  will  influence  the  material’s  properties. 

To  assist  the  modeling  effort  at  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC),  the  Center  for  Advance  Research  and 
Material  Characterization  (CARMA),  ERDC,  is  being  used  to  support  the 
polycrystalline  center-directed  research  project.  The  CARMA  is  equipped 
with  Zeiss  optical  microscopy  (OM)  workstations  and  an  environmental 
scanning  electron  microscope  (eSEM)  manufactured  by  FEI  Company, 
Hillsboro,  OR,  which  allows  for  micro-  and  nano-scale  material 
characterization.  Mechanical  property  information  at  submicron  length 
scale  is  obtained  by  a  nanoindenter  (Agilent,  Santa  Clara,  CA)  and  an  in 
situ  SEM  picoindenter  (Hysitron,  Minneapolis,  MN).  At  higher  length 
scales,  universal  testing  machines  (Instron,  Norwood,  MA)  and  load 
frames  by  MTS,  Eden  Prairie,  MN,  with  loading  capabilities  ranging  from 
5  kN  to  200  kN  provide  stress  state  and  temperature  effects  analysis  for 
quasi-static  and  cyclic  loading  conditions.  The  CARMA  is  also  capable  of 
characterizing  elemental  composition  of  materials  by  implementing  either 
a  wavelength-dispersive  X-ray  fluorescence  (XRF)  spectrometer 
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(Panalytical,  Almelo,  The  Netherlands)  or  an  energy-dispersive  X-ray 
(EDX)  spectrometer  detector  (Bruuker,  Billerica  MA)  attached  to  the 
eSEM.  The  arrangement  of  atoms  within  a  crystal  structure  may  be 
characterized  with  a  Panalytical  X-ray  diffractometer  (XRD).  Besides  the 
characterization  capabilities  at  CARMA,  collaborations  between  various 
universities  and  government  laboratories  have  been  initiated.  Carnegie- 
Mellon  University  will  assist  in  3-D  microstructure  analysis;  the  Center  for 
Advanced  Vehicular  Materials  (CAVs)  at  Mississippi  State  University 
provided  their  image  analysis  program;  the  Beckmann  Institute  at  the 
University  of  Illinois  Urbana-Champaign  provides  additional  material 
testing  and  characterization  capabilities;  and  the  Center  for  Advanced 
Microstructures  and  Devices  (CAMD)  at  Louisiana  State  University  will 
provide  in  situ  XRD  mechanical  testing  and  computed  tomography. 

This  chapter  discusses  characterization  and  testing  performed  to  date  of 
fracture  toughness,  elastic  modulus,  hardness,  grain  size,  and  crystalline 
phases  examined  for  a  commercially  available  SiC,  which  will  aid  in 
providing  the  necessary  structure-property  relationships  of  various 
polycrystalline  materials  to  aid  in  microstructure-based  model 
development,  calibration,  and  validation  efforts  at  the  ERDC. 

6.1  Material  preparation 

Laboratory  capabilities  were  demonstrated  through  tests  of  as-sintered  SiC 
commercially  available  from  Saint-Gobain  Ceramics,  Amherst,  NY, 
obtained  from  the  Army  Research  Laboratory  (ARL),  Adelphi,  MD.  The 
original  size  of  the  as-received  tile  was  4.22  by  4.22  by  0.36  in.  The 
fracture  toughness,  nanoindentation,  and  XRD  specimens  were  all 
machined  from  the  same  tile.  In  addition,  material  from  the  tile  was 
subjected  to  grain  size  analysis.  Fracture  toughness  specimens  were 
machined  using  a  Secotom-io  high-precision  cut-off  saw  with  a  diamond 
wheel  oil-cooled  blade  (Struers,  Westlake,  OH).  Machined  specimens  had 
dimensions  of  0.2  by  0.3  by  3.0  in.  with  cutting  occurring  down  the  length 
of  the  0.2-in.  sides.  Surface  grinding  occurred  on  the  0.3-in.  side  with  a 
36-micron  diamond  wheel  to  remove  surface  defects  from  manufacturing. 

After  initial  machining,  a  digital  caliper  (Mitutoyo,  Chicago,  IL)  measured 
the  dimensions  of  the  length,  width,  and  thickness  of  the  fracture 
toughness  specimen.  Three  measurements  for  each  direction  were 
obtained.  Additionally,  for  density  determinations  of  the  fracture 


ERDC  TR-12-9 


70 


toughness  specimens,  a  laboratory  balance  (Sartorius,  Bradford,  MA) 
measured  the  wet  and  dry  masses  of  the  specimens. 

Once  the  mass  and  dimensional  measurements  were  completed,  the 
fracture  toughness  specimens  were  notched  by  precision  cuts  in  the  center 
of  the  specimens  using  a  0.015-in. -thick  diamond  wheel  resulting  in  a  ratio 
of  average  notch  depth  to  specimen  depth,  a0/d,  of  0.378.  An  optical 
microscope  (Imager,  Zeiss,  Thornwood,  NY)  was  used  to  examine  the 
notch  depths  after  machining,  as  depicted  in  Figure  34. 


Noticeable  machining  damage  is  evident  for  the  specimen  in  Figure  34, 
which  is  representative  of  the  other  specimens  as  well.  As  Srinivasan  and 
Seshadri  (1980)  previously  reported,  the  cracked  grains  act  as  natural 
flaws  from  which  fractures  initiate,  obviating  the  need  for  precracking  the 
specimens  prior  to  testing. 
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The  bulk  density  of  the  fracture  toughness  specimens  was  performed  by 
the  Archimedes  water  immersion  method  according  to  ASTM  B311- 
o8(American  Society  for  Testing  and  Materials  2008).  Srinivasan  and 
Seshadri  (1980)  used  this  technique  previously  to  quantify  the  bulk 
density  of  a  sintered  SiC. 

Figure  35  represents  the  4-point  bend-loading  configuration  used  to  test 
ten  single-edge  notched-beam  (SENB)  specimens.  An  Instron  8852 
(Norwood,  MA)  load  frame  loaded  the  specimens  to  fracture  at  a 
crosshead  velocity  of  0.02  in. /min  with  load  and  displacement  data 
recorded  until  fracture. 


P  =  load 

ao  =  notch  length 
b  =  specimen  width 
d  =  specimen  height 

Li  =  loading  span 
length 

L2  =  support  span 
length 

P=load 

ao=machined  notch 
length 


P/2 


P/2 


Figure  35.  Loading  configuration  for  SENB  specimens. 


Table  2  details  the  dimensions  for  the  testing  configuration  shown  in  Fig¬ 
ure  35. 


Table  2.  Loading  and  support  span 
lengths  for  the  SENB  specimens. 


Distance  (in.) 

Loading  Span  Li 

0.6 

Support  Span  L2 

1.8 

An  FEI  Nova  eSEM  630  field  emission  gun  scanning  electron  microscope 
(FEG-SEM)  provided  fracture  surface  analysis  of  the  fracture  toughness 
specimens  while  chemical  composition  analysis  used  the  energy-dispersive 
X-ray  (EDX)  detector  attached  to  the  eSEM. 
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The  nanoindentation  specimen  was  machined  using  a  Struers  Secotom-io 
high-precision  cut-off  saw  with  diamond  wheel  oil-cooled  blade.  The 
nanoindentation  specimen  dimensions  were  0.5  by  0.5  by  0.36  in.  After 
machining,  the  specimen  was  mounted  in  1.25-in.  diameter  molds  using  a 
commercially  available  epoxy  (EpoMet,  Buehler,  Lake  Bluff,  IL).  The 
procedure  used  for  performing  stepwise  polishing  of  the  specimens  is 
outlined  in  Table  3. 


Table  3.  Polishing  of  SiC  for  Nanoindentation. 


Step 

Abrasive/Size 

Base  Speed  (rpm) 

Load  (Ibf) 

Time  (min:sec) 

1 

35-micron  diamond  pad 

240 

6 

Until  plane 

2 

9-micron  diamond  paste 

240 

6 

8:00 

3 

3-micron  diamond  paste 

120 

6 

6:00 

4 

0.05micron  colloidal  silica 

120 

6 

4:00 

The  specimens  were  ultrasonically  cleaned  with  ethanol  between  polishing 
steps  while  still  held  in  the  automatic  polishing  specimen  holder  to  ensure 
that  polishing  occurred  only  in  the  one  direction.  Once  the  polished 
surface  was  analyzed  in  the  SEM,  the  specimen  was  placed  in  a  Nano 
Indenter  G200  (Agilent,  Englewood,  CO)  with  a  Berkovich  indenter  tip. 
The  nanoindenter  performed  load  control  experiments  on  the  SiC  samples 
at  maximum  loads  of  5, 10,  50,  too,  200,  and  400  mN  using  the 
traditional  quasi-static  (QS)  method  and  the  continuous  stiffness 
measurement  (CSM)  method.  To  compare  the  results  of  these  experiments 
with  previous  work  performed  by  Guicciardi  et  al.  (2007),  the  DSI 
experiments  were  performed  at  a  loading  and  unloading  rate  of  2.5  mN/s 
with  no  hold  time.  The  CSM  experiments  were  performed  at  a  surface 
approach  velocity  of  10  nm/s,  a  strain  rate  of  0.05/s,  and  a  target 
frequency  of  45  Hz.  Both  methods  performed  a  minimum  of  ten  indents 
placed  50  um  apart.  Machine  compliance  and  thermal  drift  are 
automatically  accounted  for  with  the  G200  testing  software,  which  are 
based  on  the  Oliver  and  Pharr  (1992)  analysis  of  nanoindentation  data.  A 
standard  fused  silica  reference  sample  was  used  for  tip  area  calibration. 

The  bright  field  image  feature  of  a  Zeiss  Imager  Zim  optical  microscope  in 
reflected  mode  presented  the  images  necessary  for  grain  size  analysis. 
Specimens  were  polished  according  to  the  procedure  detailed  previously  in 
Table  3.  After  polishing,  the  sample  was  etched  for  grain  size  analysis 
using  a  modified  Murakami’s  solution  outlined  in  Table  4. 
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Table  4.  Modified  Murakami’s  solution  used  for  etching. 


Constituent 

Amount 

Distilled  Water 

100  mL 

Sodium  Hydroxide 

3g 

Potassium  Ferricyanide 

30  g 

An  optimal  etching  time  of  30  minutes  in  the  boiling  solution  yielded 
visible  grain  boundaries.  The  image  analysis  process  used  a  MATLAB 
program  developed  at  CAVs  to  post-process  the  optical  microscopy  images 
and  provided  the  statistical  data  of  the  microstructure  as  described  in 
Allison  et  al.  (2011). 

6.2  Results  and  discussion 

Calculations  of  fracture  toughness  values  implemented  the  equation  by 
Brown  and  Srawley  (1966),  which  Srinivasan  and  Seshadri  (1980)  used  for 
calculations  of  their  values  as  reported  in  Table  5.  Fracture  toughness 
values  reported  by  the  manufacturer  for  SENB  experiments  are  also 
included  in  Table  5,  although  the  equation  and  initial  notch  depth  ratio 
were  not  provided. 


Table  5.  Mechanical  property  results  for  the  fracture  toughness  specimens. 


Mechanical 

Property 

Experimentally 
Determined  Value 
(this  study) 

Manufacturer  Reported 
Value  (Saint-Gobain 
Advanced  Ceramics 

2010) 

Value  Calculated  by 
Srinivasan  and 

Seshadri  (1980) 

Density  (g/cm3) 

3.13+0.05 

3.13 

3.125-3.155 

Fracture 
toughness 
(MPa  x  m0-5) 

4.65+0.483 

4.60 

4.783+0.277 

Initial  Notch 
Depth  Ratio 

( ao/d ) 

0.378±0.003 

N/A 

0.408+0.005 

Table  5  reports  the  average  and  standard  deviation  of  10  specimens  for 
this  study.  Slight  differences  in  specimen  preparation,  geometry,  and 
material  are  likely  causes  for  the  difference  in  the  average  fracture 
toughness  values  reported.  However,  the  data  for  this  study  agree  with 
manufacturer’s  reported  values. 

An  FEI  Nova  SEM  630  FEG-SEM  provided  fracture  surface  analysis  imag¬ 
es.  Under  quasi-static  loading  the  SENB  specimen  exhibited  predominate¬ 
ly  transgranular  fracture  as  Figure  36  exhibits. 
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The  predominant  transgranular  fracture  observed  on  the  fracture  surfaces 
in  this  study  was  very  similar  to  the  transgranular  fracture  depicted  in 
Shih  et  al.  (1998)  of  a  conventional  SiC  loaded  quasi-statically  in  flexure 
that  used  boron  and  carbon  as  sintering  aids. 
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Figure  36.  Fracture  surface  of  SENB  test  specimen. 


The  indentation  modulus  of  the  SiC  sample  determined  by  the  QS  and 
CSM  methods  are  plotted  in  Figure  37  in  which  Young’s  modulus  is  plotted 
versus  the  maximum  indentation  load  applied  during  testing.  The 
modulus  was  determined  by  the  nanoindenter  TestWorks  V4.8  software, 
which  automatically  corrects  for  thermal  drift  and  machine  compliance. 

In  Figure  37,  most  noticeable  is  the  reduction  in  modulus  for  the  QS 
method  data  when  loads  exceeded  10  mN.  This  reduction  indicates  more 
of  a  bulk  modulus  response  influenced  by  microstructural  features  such  as 
porosity  instead  of  a  modulus  response  of  the  individual  phases  at  lower 
loads.  The  CSM  method  did  not  depict  this  trend;  only  the  scatter  of  the 
data  appeared  reduced  at  loads  greater  than  10  mN.  Both  methods  did 
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Figure  37.  Plot  of  Young’s  modulus  calculated  using  Oliver  and  Pharr  (1992)  versus  peak 
indentation  load.  Bars  represent  maximum  and  minimum  values  of  each  data  set. 


Figure  38.  Load  versus  displacement  curves  of  QS  nanoindentation  method  experiments  for  a 

cracked  and  uncracked  test. 
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depict  larger  amounts  of  scatter  at  maximum  loads  below  10  mN  with  the 
exception  of  the  CSM  data  at  400  mN.  The  CSM  data  appears  to  show 
much  more  stability  than  the  QS  data  in  the  average  data  points  for  the 
material  across  the  various  maximum  loads.  The  reported  bulk  modulus  of 
elasticity  from  the  manufacturer’s  data  sheet  is  410  GPa.  The  manufactur¬ 
er  did  not  report  how  the  modulus  value  was  measured.  Both  QS  and  CSM 
method  averages  appear  to  overpredict  the  modulus  reported  by  the  man¬ 
ufacturer  below  10-mN  loads;  however,  the  QS  method  did  capture  the 
manufacturer’s  reported  value  in  the  scatter  of  the  data  at  5  mN  while  both 
captured  the  reported  value  at  10  mN.  This  trend  suggests  that  at  lower 
loads  the  methods  were  not  capturing  a  bulk  modulus  of  elasticity  value. 

At  loads  above  10  mN  the  scatter  in  the  data  decreased  for  both  methods. 
At  50-  and  100-mN  loads,  only  the  QS  method  captured  the  manufactur¬ 
er’s  reported  bulk  modulus  while  the  data  for  the  CSM  method  was  slightly 
increased.  For  maximum  loads  of  200  and  400  mN  the  QS  data  exhibited 
lower  modulus  values  than  the  CSM  method  data.  The  deviation  between 
the  CSM  and  QS  data  at  higher  loads  might  be  due  to  the  fact  that  the  QS 
method  provides  data  in  observable  steps,  allowing  experiments  with 
cracking  to  be  identified  and  removed  from  the  data  sets  prior  to  making 
statistical  comparison.  Figure  38  depicts  load-unload  curves  for  the  QS 
method  with  cracked  and  uncracked  responses. 

A  noticeable  step  in  the  load  displacement  curve  is  visible  for  the  cracked 
test  in  Figure  38.  Analyzing  every  load  displacement  curve  for  the  QS 
method  data  allows  researchers  to  determine  invalid  tests  that  do  not 
represent  valid  measurement  of  modulus  for  intact  material. 

Hardness  results  are  reported  for  the  QS  and  CSM  methods  in  Figure  39. 
The  hardness  was  determined  by  the  nanoindenter  TestWorks  V4.8 
software,  which  automatically  corrects  for  thermal  drift  and  machine 
compliance. 

As  observed  with  the  modulus  data  reported  in  Figure  37,  the  hardness 
data  shown  in  Figure  39  exhibits  less  scatter  at  loads  above  10  mN.  Unlike 
the  modulus  data,  the  hardness  data  exhibit  much  greater  overlap  between 
values  obtained  by  the  QS  and  CSM  methods.  In  addition,  the  average 
CSM  results  display  less  dependence  on  load  than  the  average  QS  results. 

Figure  40  provides  a  post-processed  image  of  the  etched  SiC  grains  after 
using  the  MSU  image  analyzer  program.  The  grains  are  shown  as  gray  ob- 
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jects  in  the  image,  while  the  black  objects  are  pores.  Only  the  gray  objects 
are  analyzed  for  grain  size  measurements.  A  total  of  five  images  were  ana¬ 
lyzed  for  grain  size  measurements. 

The  average  grain  size  diameter  determined  from  the  image  analysis  soft¬ 
ware  was  6.5  pm,  which  correlates  to  the  manufactured  reported  grain  size 
of  4-10  pm. 
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Figure  39.  Oliver  and  Pharr  (1992)  calculated  hardness  versus  peak  indentation  load.  Bars 
represent  maximum  and  minimum  values  of  each  data  set. 
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Figure  40.  XRD  pattern  of  commercially  available  SiC. 
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posite  that  has  fracture  toughness  and  tensile  strength  approximately  5  times  that  of  existing  polycrystalline  ceramics,  such  as  silicon 
carbide  or  boron  carbide.  If  such  a  ceramic  composite  were  developed,  then,  based  on  current  strength-to-weight  and  stiffness-to-weight 
ratios,  it  could  replace  steel  and  aluminum  for  most  structural  applications  with  an  attendant  two-thirds  reduction  in  weight.  This  would 
have  enormous  impact  on  Army  portable  protective  structures,  equipment,  and  logistics.  Key  to  this  development  is  the  growing  capa¬ 
bility  in  numerical  simulations  to  predict  material  behavior  based  on  atomic  and  crystalline  morphology.  Such  simulations  provide  new 
insight  into  the  causal  relationships  between  material  structure  and  material  behavior.  The  simulations  can  guide  both  polycrystalline 
material  design  and  synthesis  methods  such  as  sintering. 
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